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ABSTRACT 


The  principal  objective  of  this  research  is  to  obtain  an  improved 
understanding  of  the  structure  of  laminar  flames  burning  hydrocarbon  fuels, 
alcohol  fuels,  and  blends  of  hydrocarbon  and  alcohol  fuels.  This  work  attempts 
to  support  the  U.  S.  Army's  research  and  development  programs  in  alternative 

fuels.  Since,  there  is  convincing  evidence  that  turbulent  reacting  flows  can 
be  modelled  as  statistical  ensemble  of  laminar  flamclets,  an  experimental, 
numerical  and  analytical  study  was  undertaken  to  clarify  the  structure  and 
mechanisms  of  extinction  of  laminar  diffusion  flamclets,  laminar  premixed 
flamclets  and  laminar  partially  premixed  flamclets. 

Results  of  the  study  on  diffusion  flames  burning  hydrocarbon  and 
alcohol  fuel  show  that  the  basic  structure  of  these  flames  arc  essentially 

similar  and  can  be  modelled  by  use  of  four  overall  reactions.  Only  one  of  these 
overall  reactions  is  specific  to  the  fuel  being  considered,  while  the  other  three 
reactions  are  those  which  describe  the  oxidation  of  H2  and  CO  to  H2O  and  C02-  It 
was  found  that  chemical  reactions  occur  in  roughly  two  distinct  layers  which 
are  termed  the  inner  layer  and  the  oxidation  layer.  In  the  inner  layer  the  fuel 

is  attacked  by  the  radicals  to  form  H2  and  CO,  which  subsequently  oxidize  in 

the  oxidation  layer  to  H2O  and  CO2.  This  description  of  the  flame  structure 
considerably  simplifies  the  analysis.  It  was  also  found  that  the  structure 
premixed  flames  and  partially  premixed  flames  consist  of  these  distinct  layers. 
This  finding  implies  that  accurate  analytical  description  of  the  combustion 
process  in  la  ninar  flamclets  are  possible,  and  the  results  can  be  extended  to 
turbulent  flames. 
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CHAPTER  I 


INTRODUCTION 

The  objective  of  this  research  was  to  support  the  U.  S.  Army  research  and 
development  program  in  alternative  fuels.  The  research  was  performed  in 
collaboration  with  Professor  N.  Peters  at  the  Institut  of  Tcchnischc  Mcchanik, 
Rhcinish-Wcstfalische  Tcchnischc  Hochschule  at  Aachen,  West  Germany, 
Professor  C.  Trevino  at  Dcpto  Fluidos  y  Termica,  Facultad  dc  Ingenieria,  UNAM, 
Mexico,  and  Professor  M.  D.  Smooke  at  the  Department  of  Mechanical 
Engineering,  Yale  University,  New  Haven,  Connecticut. 

The  principal  objective  of  this  research  is  to  obtain  an  improved 
understanding  of  the  combustion  of  hydrocarbon  and  alcohol  fuels.  Simple, 
but  fundamentally  sound  methods  for  evaluating  the  combustion  of 
hydrocarbon  fuels  have  also  been  developed.  The  structure  of  premixed 
flames  and  diffusion  flames  were  considered.  In  the  previous  contract  entitled 
"Mechanisms  of  Combustion  of  Hydrocarbon/Alcohol  Fuel  Blends", 

experiments  were  performed  to  characterize  the  structure  of  diffusion  flames 
burning  hydrocarbon  fuels,  alcohol  fuels  and  blends  of  hydrocarbon  and 
alcohol  fuels.  A  major  effort  in  the  present  program  involved  interpretation 
of  these  experimental  results.  Detailed  numerical  calculations  and  asymptotic 

analysis  was  performed  to  interpret  these  experimental  results.  The 
techniques  developed  here  can  be  readily  extended  to  most  hydrocarbon  fuels 
md  fuel  blends. 

Since  most  of  the  results  of  this  research  have  ecu  published  in  the 
open  unclassified  literature,  only  the  abstract  of  these  published  work  is  given 
here,  and  the  reader  is  referred  to  the  literature  for  further  details.  Results  of 
research  which  hove  been  submitted  for  publication  arc  reproduced  here. 

In  Chapter  II  results  of  numerical  studies  on  heptane-air  diffusion  flames 
arc  reported.  It  was  found  that  the  chemical  kinetic  mechanism  for  oxidation 
of  heptane  can  be  reduced  to  four  overall  reactions.  In  Chapter  III  a 

corresponding  numerical  study  was  performed  on  methanol-air  diffusion 

flames,  and  the  results  show  that  the  overall  chemical  kinetic  mechanism  of 


oxidation  of  methanol  can  also  be  reduced  either  to  five  overall  reactions. 
Although  these  numerical  studies  provide  valuable  information  concerning 
the  structure  of  these  flames  they  require  substantial  computer  time.  In 
addition,  numerical  techniques  may  mask  the  underlying  physics  of  the 
phenomena.  Asymptotic  analysis  on  the  other  hand  require  minimal 

computer  time,  and  often  clarify  the  essential  physics  of  the  phenomena.  In 
order  to  develop  asymptotic  techniques  for  analyzing  the  structure  of 
hydrocarbon  flames,  studies  were  made  on  simple  hydrocarbon  fuels  such  as 
methane  and  the  results  arc  summarized  in  Chapter  IV.  Since  turbulent 
combustion  can  be  modelled  as  a  statistical  ensemble  of  laminar  diffusion 
flamclets,  laminar  premixed  flamclcts,  and  laminar  partially  premixed 
flamelcts,  experimental,  numerical  and  asymptotic  studies  were  performed  on 
laminar  premixed  flames  and  they  are  reported  in  Chapter  V.  Similarly  in 
Chapter  VI  studies  on  partially  premixed  flames  arc  reported.  It  was  found  that 
that  hydrogen-oxygen  reactions  play  an  in  the  combustion  of  hydrocarbon 
and  alcohol  fuels.  Therefore,  a  numerical  and  analytical  study  was  undertaken 
to  determine  the  structure  of  hydrogen-air  diffusion  flames  and  the  results 
are  summarized  in  Chapter  VII.  Since,  the  characteristic  Lewis  Number  of 
hydrogen  is  small,  differential  diffusion  of  this  species  plays  an  important 
role  in  determining  the  structure  of  the  flame.  Hence,  in  Chapter  VII,  results 
of  analytical  studies  of  the  influence  of  Lewis  Number  of  the  reactants  on  the 
asymptotic  structure  of  countcrflow  diffusion  flames  arc  also  summarized. 


CHAPTER  II  HEPTANE  -  AIR  DIFFUSION  FLAMES 


II  .1  STRUCTURE  OF  HEPTANE-AIR  COUNTERFLOW  DIFFUSION 

FLAMES 

The  Supercomputer  located  at  the  NSF  Center  at  San  Diego  was  used  to 
characterize  the  structure  of  heptane-air  diffusion  flames.  The  computer 
program  and  formulation  for  the  numerical  problem  used  in  these 
calculations  were  developed  by  Professor  M.  D.  Smookc  at  Yale  University.  The 
chemical  kinetic  mechanism  for  oxidation  of  heptane  used  in  this  paper 
consists  of  forty-two  elementary  reactions  involving  eighteen  species.  The 
balance  equations  for  mass,  momentum,  energy,  and  species  were  integrated 
numerically,  and  converged  solutions  were  typically  obtained  using 
approximately  four  hours  of  CPU  time.  The  results  of  numerical  calculations 
were  compared  with  experimental  measurements,  and  good  agreement  was 
obtained. 

The  research  summarized  above  has  been  published  in  the  Proceedings  of 
the  Fourth  International  Conference  on  Supercomputer  and  Third  World 
Supercomputer  Exhibition,  Volume  II,  Edited  by  Professor  Lana  P.  Kartashev 
and  Dr.  Steven  I.  Kartashev,  International  Supercomputing  Institute,  Inc,  Suite 
B-309,  3000-34th  Street,  South.  St.  Pctersburgh,  Florida  33711,  pp  490-495,  1989. 
The  coauthors  of  this  publication  were  Ms.  Mary  Bui-Pham  and  Professor  F.  A. 
Williams 
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D.  2  COMPARISON  BETWEEN  EXPERIMENTAL  MEASUREMENTS  AND 


NUMERICAL  CALCULATIONS  OF  THE  STRUCTURE  OF 
HEPTANE-AIR  DIFFUSION  FLAMES 

ABSTRACT 

Detailed  numerical  calculations  are  performed  to  determine  the  structure  of  heptane-air 
diffusion  flames,  and  the  results  are  compared  with  experimental  measurements.  The 
configuration  used  is  the  diffusion  flame  stabilized  in  the  vicinity  of  a  stagnation  plane, 
which  is  formed  by  directing  an  oxidizing  gas  flow  onto  the  vaporizing  surface  of  a  pool  of 
heptane.  Profiles  of  the  concentration  of  various  stable  species  and  of  the  temperature 
have  been  measured  by  gas  chromatography  and  by  thermocouples,  respectively.  To 
evaluate  the  influence  of  strain  on  the  structure  of  the  flame,  the  measurements  taken  at  a 
fixed  composition  of  the  oxidizer  stream  and  at  two  values  of  the  strain  rate  were  chosen 
for  comparison  with  the  calculated  results.  The  computations  were  performed  using  a 
chemical  kinetic  mechanism  consisting  of  forty-two  elementary  reactions  involving  eighteen 
species.  To  simplify  the  chemical  kinetic  mechanism,  it  was  assumed  that  heptane  is 
attacked  by  radicals  to  form  the  heptyl  radical  whose  decomposition  to  CH3  and  C3H6  is 
represented  by  a  one-step  overall  reaction.  Good  agreement  was  obtained  between  the 
results  of  the  detailed  numerical  calculations  and  the  experimental  measurements.  It  was 
found  that  aspects  of  the  structure  of  heptane-air  diffusion  flames  relevant  to  asymptotic 
descriptions  are  similar  to  those  of  methane-air  diffusion  flames.  Hence,  the  previously 
developed  techniques  for  asymptotic  analysis  of  the  structure  of  methane-air  flames  can  be 
used  for  analyzing  the  structure  of  heptane-air  flames. 


1.  INTRODUCTION 


Motivated  by  recent  success  in  numerical  and  analytical  modeling  of  the  structure  of 
laminar  flames  burning  simple  hydrocarbon  fuels  such  as  methane  and  propane,  Warnatz 
(1984),  Esser  et  al.  (1985),  and  Westbrook  et  al.  (1988)  have  studied  the  combustion  of 
more  complex  hydrocarbon  fuels  such  as  n-heptane  and  iso-octane  which  are  often  used  as 
reference  fuels.  These  previous  investigations  have  been  concerned  with  premixed 
systems.  In  this  paper,  results  of  numerical  calculations  of  heptane-air  diffusion  flames 
structure  are  reported  and  compared  with  previous  experimental  measurements  made  by 
Kent  and  Williams  (1974),  Seshadri  (1977),  and  Hamins  and  Seshadri  (1987).  The 
calculations  were  performed  on  the  Cray  X-MP  located  at  the  NSF  Center  at  San  Diego 
using  a  program  developed  by  Smooke  (1982). 

To  use  a  nearly  complete  description  of  the  chemical  kinetics  of  oxidation  of  heptane 
in  the  calculations  would  be  difficult  because  it  involves  hundreds  of  elementary  reactions 
containing  numerous  species,  some  of  which  would  have  several  isomeric  structures.  For 
example,  Warnatz  (1984)  suggested  that  there  may  exist  39  isomeric  structures  for  the 
principal  alkyl  radical  which  is  formed  from  heptane.  Since  the  structures  of  steady 
diffusion  flames  are  of  interest  here,  only  the  high  temperature  oxidation  of  heptane  is 
relevant.  Warnatz  (1984)  has  also  outlined  a  procedure  for  characterizing  the  chemical 
kinetic  mechanism  of  oxidation  of  higher  alkanes  including  heptane.  It  is  presumed  that 
alkanes  are  attacked  by  the  radicals  H,  O,  and  OH  to  form  the  alkyl  radical  which 
decomposes  to  smaller  radicals  by  fast  thermal  elimination  of  alkenes  to  eventually  form 
CH3  and  C3H6;  reactions  leading  to  these  compounds  are  not  presumed  to  be  rate  limiting. 
This  approximation  considerably  simplifies  the  chemistry  because  the  oxidation  of  CH3  and 
C3H6  is  better  understood.  Using  this  approximation,  Warnatz  (1984)  calculated  the 
burning  velocities  of  premixed  flames  for  a  number  of  higher  hydrocarbon  fuels  including 


heptane  and  found  good  agreement  between  the  calculated  values  and  experimental 
measurements.  The  calculated  burning  velocities  were  found  to  be  relatively  insensitive  to 
variations,  within  reasonable  limits,  of  the  rate  of  decomposition  of  the  principal  alkyl 
radical  to  CH3  and  C3H6.  Therefore,  in  these  calculations,  the  decomposition  of  C7H15 
was  represented  by  a  one-step  reaction  to  form  CH3  and  CsHg,  and  the  rate  of  this  reaction 
was  chosen  to  be  equal  to  that  of  decomposition  of  the  n-butyl  radical.  Studies  of  this  type 
are  useful  in  identifying  key  reactions  characterizing  the  structure  and  mechanisms  of 
oxidation  of  hydrocarbon  fuels.  The  results  obtained  here  are  useful  for  future  asymptotic 
analyses  using  multi-step  overall  chemical  kinetic  models  similar  to  those  performed  by 
Seshadri  and  Peters  (1988)  for  methane-air  diffusion  flames. 

2.  DESCRIPTION  OF  THE  EXPERIMENTAL  MEASUREMENTS  AND 
FORMULATION  OF  THE  NUMERICAL  PROBLEM 

2.1  Description  of  the  Experimental  Measurements: 

Experiments  were  performed  by  Kent  and  Williams  (1974),  Seshadri  (1977),  and 
Hamins  and  Seshadri  (1987)  to  determine  the  structure  of  diffusion  flames  stabilized  in  the 
stagnation  point  boundary  layer  by  directing  an  oxidizing  gas  stream  vertically  downward 
onto  the  vaporizing  surface  of  pools  of  heptane.  The  oxidizer  duct  was  located  1  cm  above 
the  fuel  surface;  a  diffusion  flame  can  be  indefinitely  stabilized  in  this  configuration. 
Seshadri  and  Williams  (1978)  suggested  that  the  strain  experienced  by  the  flame  sheet  can 
be  approximated  if  the  flow  velocity  at  the  exit  of  the  oxidizer  duct  is  known,  and  for  given 
values  of  strain  and  composition  of  the  oxidizing  gas,  the  fuel  burning  rate  and  the  surface 
temperature  can  be  calculated  if  the  thermodynamics  of  vaporization  and  oxidation  are 
known. 
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Concentration  profiles  of  stable  species  were  measured  using  gas  chromatographic 
techniques,  and  temperature  profiles  were  measured  by  thermocouples.  The  two  sets  of 
experimental  measurements  chosen  for  comparison  with  numerical  calculations  are:  1)  Y(*. 
=  0.185,  U=  30.5  cm/s,  and  2)  Yq20O  =  0.185,  U=  83.7  cm/s,  where  Y0^  and  U  are  the 
mass  fraction  and  axial  velocity  of  the  oxidizing  gas  stream  at  the  exit  of  the  duct, 
respectively.  Flames  at  two  different  values  of  strain  with  the  same  values  of  are 
considered  here  to  clarify  the  influence  of  strain  on  the  structure  of  the  flame.  Since  the 
theoretically  calculated  strain  rate  is  directly  proportional  to  the  velocity  of  the  oxidizing  gas 
at  the  exit  of  the  duct,  the  increase  in  the  strain  rate  can  be  obtained  by  increasing  the 
oxidizer  stream  velocity.  Experimental  data  for  U  =  30.5  cm/s  were  collected  by  Kent  and 
Williams  (1974),  and  those  for  U=  83.7  cm/s  are  shown  in  Figs.  8  and  9.  Further  details 
of  the  experimental  apparatus  and  procedure  for  making  these  measurements  are  outlined 
by  Kent  and  Williams  (1974),  Seshadri  (1977),  and  Hamins  and  Seshadri  (1987). 

2.2  Chemical  Kinetic  Mechanism: 

Table  1  shows  the  chemical  kinetic  mechanism  used  in  the  calculations.  The  rate 
constants  are  assumed  to  be  in  the  Arrhenius  form,  kj  =  Aj  Tnj  exp  (-Ej  /  (RT)),  where  Aj 
is  the  frequency  factor,  nj  is  the  temperature  exponent,  and  Ej  is  the  activation  energy.  The 
data  shown  in  Table  1  represent  the  forward  rate  of  the  reactions;  the  backward  rates  are 
calculated  from  the  equilibrium  constant. 

The  data  selected  for  the  rate  constants  for  reactions  1-33  were  obtained  from 
published  work  of  Wamatz  (1984),  Westbrook  and  Dryer  (1984),  and  Miller  et  al  (1984), 
and  they  are  similar  to  those  used  previously  by  Smooke  et  al.  (1986),  Puri  et  al.  (1987), 
and  Smooke  et  al.  (1988)  for  methane-air  flames  and  partially  premixed  flames.  The 
results  of  the  numerical  calculations  using  these  rates  were  found  to  agree  well  with 
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experimental  measurements.  The  data  selected  for  the  rates  of  reactions  34-96  are  those 
recommended  by  Wamatz  (1984). 

2.3  Formulation  of  the  Numerical  Problem: 

The  formulation  of  the  numerical  problem  is  described  in  detail  by  Smooke  et  al. 
(1986);  hence,  only  the  main  features  are  outlined  here.  Let  r  and  z  denote  the  independent 
spatial  coordinates  in  the  radial  and  axial  directions,  respectively.  The  primary  oxidizer 
stream  is  located  at  z  =  LM,  and  the  surface  of  the  liquid  fuel  at  z  =  0.  The  experimental 
flow-field  consists  of  an  inner,  viscous  flow  region  near  the  fuel  surface  which  extends 
slightly  beyond  the  stagnation  plane  into  the  oxidizer  side  and  an  outer,  inviscid,  and 
rotational  region  extending  from  the  oxidizer  side  of  the  stagnation  plane  to  the  exit  of  the 
oxidizer  duct,  where  the  tangential  component  of  the  flow  velocity  vanishes.  Hence,  in  the 
outer  flow,  the  strain  rate  a  is  a  function  of  the  axial  coordinate  z,  as  suggested  by  Seshadri 
and  Williams  (1978).  However,  in  the  computational  model  used  here,  the  outer  flow  is 
assumed  to  be  inviscid  and  irrotational,  and  the  tangential  component  of  the  flow  velocity 
uM  at  the  oxidizer  duct  is  not  zero.  Thus,  in  the  outer  flow,  the  value  of  a  is  constant,  and 
the  tangential  and  normal  component,  u  and  v,  at  the  edge  of  the  boundary  layer  can  be 
written  as  uM  =  ar,  and  v«.  =  -2az,  where  <*>  denotes  conditions  at  the  primary  oxidizer 
stream  at  z  =  LM.  Hence,  the  computational  model  used  for  the  outer  flow  is  different  than 
that  encountered  in  the  experiment  and  results  in  a  slight  shift  in  the  flame  location.  A  more 
accurate  description  of  the  flow  field  incorporating  the  rotational  aspect  of  the  outer  flow  is 
described  by  Kee  et  al.  (1988). 

The  notations  f  =  u/uoo  and  M  =  pv  are  introduced,  where  f  is  related  to  the 
derivative  of  a  modified  stream  function,  as  suggested  by  Dixon-Lewis  et  al.  (1984),  and  p 
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is  the  gas  density.  The  governing  boundary  layer  equations  for  mass,  momentum,  chemical 
species  and  energy  in  cylindrical  coordinates  are: 


4-M.  +  2apf=  0 

+a(p~-pff|2)=0 

--MpYkVk)-M^  +  WkWk  =  0,  k  =  1,  2, K 
dz  dz 


-CpMdX-  ipYkVkcpkil-  £ 

dz  k=1  dz  k=, 


wkWkhk  =  0 


and  the  system  is  closed  with  the  ideal  gas  law 

p=pW/  (RT) 

In  the  above  set  of  equations,  T  denotes  the  temperature;  Yk,  the  mass  fraction  of  the 
k™  species;  p,  the  pressure;  wk,  the  molar  rate  of  production  of  the  species;  Wk,  the 
molecular  weight  of  the  k™  species;  W,  the  mean  molecular  weight  of  the  mixture;  R,  the 
universal  gas  constant;  X,  the  thermal  conductivity  of  the  mixture;  Cp,  the  heat  capacity  of 
the  mixture;  Cpk,  the  heat  capacity  of  the  k^  species;  hk,  the  specific  enthalpy  of  the  k1^ 
species;  ji,  the  viscosity  of  the  mixture,  and  Vk  is  the  diffusion  velocity  of  the  k1*1  species 
in  the  z  direction.  Equations  for  determining  Vk,  and  wk  are  outlined  by  Kee  et  al.  (1983). 
Since  the  mass  flux  and  composition  of  the  oxidizing  gas  at  the  exit  of  the  duct  are  known, 
the  undetermined  strain  rate  can  be  calculated  as  an  eigenvalue  of  the  problem  by 
introducing  the  differential  equation  da/dz  =  0. 


To  complete  the  specifications  of  the  problem,  appropriate  boundary  conditions  must 
be  imposed.  At  the  oxidizer  jet  located  at  z  =  L^,  the  conditions  consistent  with  the  model 
used  for  the  outer  flow  are: 

f  =  1  ;  Y02  =  Y02oo,  Yk  =  0,  k*02,N2;M  =  M.;T  =  T.. 

At  the  surface  of  the  liquid  fuel,  the  tangential  component  of  the  flow  velocity  is 
presumed  to  be  zero  (no  slip),  and  the  appropriate  interface  balance  conditions  are: 

f*  w  =  0 

PwYkwVkw  +  MwYkw  =  0,  k  =  1,  2,...,K;  k  ^  F 

PwYfwVfw  +  MwYpw  =  m=  Mw 


Here  m  is  the  unknown  mass  burning  rate  and  will  be  determined  as  part  zi  the  solution;  L 
is  the  latent  heat  of  vaporization  of  the  fuel ;  F  refers  to  the  fuel,  and  w  refers  to  ?cndiMor.$ 
on  the  gas  side  of  the  liquid  gas  interface.  To  determine  the  surface  teotpcrabae  of  the  fuel, 
it  is  presumed  that  there  exists  gas-vapor  equilibrium  at  the  interface,  which  is  an  •.  •  ileiu 
approximation  for  the  values  of  m  encountered  here. 

The  governing  equations  form  a  system  of  nonlinear  two-point  boundary  value 
problems.  The  solution  procedure  as  developed  by  Smooke  (1985)  employs  a  flame  sheet 
starting  model  with  an  adaptive  time-integration  and  finite  difference  method.  Initially,  the 
surface  temperature  of  the  fuel  Tw  was  assumed  to  be  equal  to  the  normal  boiling  point, 
and  a  converged  solution  was  obtained  assuming  a  reasonable  value  for  the  burning  rate  m 
using  the  boundary  conditions  at  the  oxidizer  duct  exit  and  all  of  the  liquid-gas  interface 
balance  equations,  except  that  for  energy  balance.  To  calculate  Tw  for  vapor-liquid 


equilibrium  at  the  interface,  the  computed  value  of  the  concentration  of  the  fuel  in  the 
gaseous  form  Yfw  and  the  thermodynamic  data  collected  by  Rossini  et  al.  (1947)  were 
used.  A  revised  value  for  rft  was  calculated  using  the  energy  balance  condition,  the 
computed  values  of  the  temperature  gradient  (dT/dz)w,  and  X  on  the  gas  side  of  the 
interface.  Calculations  were  repeated  using  these  revised  values  of  Tw  and  ih  until  the 
changes  in  their  values  were  less  than  2%. 

3.  RESULTS  AND  DISCUSSION 

3.1  Simplification  of  the  Chemical  Kinetic  Mechanism: 

To  simplify  the  chemical  kinetic  mechanism,  the  numerical  calculations  were  initially 
performed  to  determine  the  structure  of  diffusion  flames  stabilized  near  the  stagnation 
region  formed  between  two  steady,  laminar,  infinitely  wide,  axisymmetric  counterflowing 
jets  of  heptane  vapor  diluted  with  nitrogen  and  air.  The  calculations  were  performed  using 
the  chemical  kinetic  mechanism  and  rate  data  shown  in  Table  1  and  for  a  value  of  the  strain 
rate,  a,  equal  to  9.84  s'1.  Results  of  the  computations  are  interpreted  using  a  suitably 
defined  mixture  fraction,  Z,  as  the  independent  variable.  This  permits  analysis  to  be 
performed  without  reference  to  specific  flow  configurations.  There  exist  several  definitions 
of  the  mixture  fraction,  that  introduced  by  Peters  (1984)  which  is  based  on  the  elements 
present  in  the  fuel,  that  introduced  by  Bilger  (1984)  which  is  based  on  the  elements  carbon 
and  hydrogen  present  in  the  fuel  and  oxygen  in  the  oxidizer  stream,  and  that  introduced  by 
Tangirala  et  al.  (1989)  which  is  based  on  an  inert  species  in  the  system  such  as  nitrogen. 
Although  the  definition  introduced  by  Bilger  (1984)  is  employed  here,  the  qualitative  and 
quantitative  results  obtained  using  other  definitions  were  essentially  the  same  for  heptane- 
air  flames.  The  mixture  fraction  Z  can  be  expressed  by  the  relation 


exactly  equal  to  zero  at  the  injection  plane;  thus,  the  value  of  a  can  be  expected  to  lie 
between  24.8  and  30.5  s'1,  although  probably  closer  to  the  latter  value.  Hence,  the  slight 
shift  between  the  computed  and  measured  peak  temperature  may  be  at  least  partially 
attributed  to  inaccuracies  in  the  flow  field  descriptic  Experimental  uncertainty  in  probe 
positioning  may  also  play  a  role. 

Asymptotic  analysis  of  the  structure  of  diffusion  flames,  as  discussed  by  Williams 
(1985)  and  Rossini  (1947),  shows  that  the  outer  flow  affects  the  flame  location,  and  the 
chemistry  is  confined  to  the  inner  zone  where  there  exists  a  diffusive-reactive  balance. 
Hence,  if  the  results  are  plotted  using  the  defined  mixture  fraction  Z  as  the  independent 
variable,  then  better  agreement  can  be  expected. 

In  Fig.  7,  the  results  shown  in  Fig.  6  are  replotted  versus  Z  and  show  that  the 
measured  and  calculated  peak  flame  temperature  occur  at  nearly  the  same  value  of  Z  = 
0.055.  Similarly,  the  measured  and  calculated  profiles  of  C7H16,  O2,  and  H2  are  aligned 
and  show  good  agreement.  However,  the  computed  values  of  CO2,  H2O,  and  CO  are 
higher  than  the  corresponding  measured  values,  as  is  also  seen  in  Fig.  6.  Since  the 
reaction  rates  along  the  path  where  H2  and  CO  are  oxidized  to  H2O  and  CO2  are  reasonably 
well  understood,  a  potential  source  of  the  discrepancies  could  be  inaccuracies  in  the 
reaction  rates  along  the  paths,  shown  in  Fig.  4,  where  C7H16  is  converted  to  CH3;  a  lower 
rate  along  the  latter  path  could  provide  better  agreement.  When  the  reaction  rates  along  this 
path  were  examined,  we  found  that  the  rate  of  reaction  37,  C3H6  +  O  -->  2  CH3  +  CO,  is 
slowest  among  the  fastest  rates.  Therefore,  calculations  were  made  using  a  frequency 
factor  for  reaction  37  which  is  10  times  lower  than  that  originally  used.  However,  no 
significant  changes  in  peak  concentration  for  the  stable  species  were  detected,  which 
indicates  that  the  rates  for  these  early  steps  are  not  responsible  for  the  discrepancies.  As  a 
further  exploratory  check,  numerical  calculations  were  repeated  with  the  frequency  factor  of 
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here  aim  is  the  stoichiometric  coefficient  denoting  the  number  of  atoms  of  element  m  in  a 
molecule  of  species  i;  Ym  p  is  the  value  of  Ym  in  the  fuel,  and  Yo.o  is  the  value  of  Yo  in 
the  ambient  oxidizer  stream. 

The  dashed  lines  in  Figs.  1  and  2  show  profiles  of  temperature  and  concentration  of 
the  major  stable  species  (C7H16,  O2,  CO2,  H2O,  H2  and  CO)  calculated  using  all  the  96 
elementary  chemical  reactions  shown  in  Table  1.  The  dashed  lines  in  Fig.  3  show  similar 
profiles  for  selected  radicals  (H,  OH,  O  and  HO2).  The  profiles  of  the  reaction  rates  of  the 
various  elementary  reactions  were  also  calculated.  The  results  of  these  calculations  were 
used  to  deduce  the  principal  path  of  oxidation  of  heptane  which  is  shown  schematically  in 
Fig.  4.  The  fuel  is  attacked  by  the  radicals  H,  OH,  and  O  to  form  C7H15  which 
decomposes  to  form  CH3  and  C3H6.  Propene  is  then  attacked  by  radicals  to  form  CH3  and 
CH3HCO,  and  CH3HCO  reacts  with  radicals  to  form  additional  CH3.  Subsequent 
oxidation  proceeds  similarly  to  that  in  methane-air  diffusion  flames.  CH3  reacts  with 
radicals  to  form  formaldehyde  which  produces  HCO  followed  by  the  formation  of  H2  and 
CO,  which  are  oxidized  to  form  H2O  and  CO2. 

To  simplify  the  chemical  kinetic  mechanism,  all  species  which  are  not  formed  or 
destroyed  along  the  principal  path  of  oxidation  of  heptane  shown  in  Fig.  4  were  removed. 
Consequently,  reactions  43-96  which  involve  the  species  CH4,  CH2,  CH,  C2H,  C2H6, 
C2H5,  C2H4,  C2H3,  C2H2,  C3H8,  C3H7,  CH2CO,  and  C2HO  were  eliminated.  The  solid 


lines  in  Figs.  1,  2,  and  3  represent  profiles  calculated  using  reaction  1-42  shown  in  Table 
1,  and  they  agree  reasonably  well  with  the  profiles  calculated  using  the  complete  set  of 
reactions  1-96.  Hence,  further  calculations  were  performed  using  reactions  1-42. 

3.2  Comparison  between  Numerical  Calculations  and  Experimental  Measurements: 

Numerical  calculations  were  performed  on  a  domain  of  1.00  cm.  For  the  cases  of  U 
=  30.5  cm/s  and  U  =  83.7  cm/s,  the  mass  fluxes  at  the  exit  of  the  duct,  M*.  were  0.0358 
gm/cm2/s  and  0.098  gm/cm2/s,  respectively,  and  the  temperature  of  the  ambient  oxidizing 
gas  stream  was  =  298  K. 

Fig,  5  shows  the  calculated  profiles  of  the  axial  velocity  along  the  axis  of  symmetry. 
In  both  cases,  the  stagnation  plane  lies  close  to  the  fuel  surface  and  moves  closer  with 
increasing  strain.  In  Fig.  6,  experimentally  measured  profiles  (shown  as  points)  for 
temperature  and  the  stable  species  C7H16,  O2,  CO2,  H2O,  H2,  and  CO  are  compared  with 
the  calculated  values  (shown  as  lines)  for  U  =  30.5  cm/s,  and  the  shapes  show  good 
agreement.  However,  the  computed  peak  temperature  occurs  0.4  mm  nearer  the  oxidizei 
side  than  the  measured  peak  temperature;  a  similar  shift  is  also  observed  for  the  species 
profiles.  A  similar  but  larger  effect  was  found  previously  for  methane-air  diffusion  flames 
by  Smooke  et  al.  (1986)  and  Puri  et  al.  (1987)  and  is  partially  attributed,  as  mentioned  in 
section  2.3,  to  the  difference  between  the  experimental  and  computational  flow  field.  The 
calculated  strain  rate  was  24.68  s"1.  However,  Seshadri  and  Williams  (1978)  showed  that, 
if  the  outer  flow  is  assumed  to  be  inviscid  and  rotational,  then  a  =  U/L  =  30.5  s_1  at  the 
stagnation  plane  which  is  presumed  to  coincide  with  the  surface  of  the  fuel.  If  the 
calculations  are  repeated  with  this  higher  value  of  a,  then  the  peak  temperature  decreases  by 
15  K  and  occurs  0.275  mm  nearer  the  fuel  side  than  the  measured  peak  temperature. 
However,  the  tangential  component  of  the  oxidizer  velocity  in  the  experiments  may  not  be 


the  rate  of  reaction  25,  HCO  +  M  -->  K  +  CO  +  M,  lowered  to  1.6  x  1013,  and  the 
calculated  peak  concentration  of  CO  shows  better  agreement;  however,  the  discrepancies 
for  CO2,  H2O  still  remain.  Hence,  it  appears  that  the  discrepancies  can  be  attributed  mainly 
to  experimental  inaccuracies,  especially  since  the  differences  are  comparable  in  magnitude 
with  estimated  experimental  uncertainties. 

Peters  (1984)  showed  from  an  asymptotic  analysis  of  the  outer  zone  that  if  the 
chemical  reaction  was  presumed  to  occur  as  a  one-step  irreversible  process,  C7H16  +  1 1  O2 
— >  7  CO2  +  8  H2O,  and  if  the  Lewis  number  of  all  species  are  presumed  to  be  nearly 
unity,  then  the  value  of  Z,  at  the  position  where  carbon,  hydrogen,  and  oxygen  are  in 
stoichiometric  proportions,  is  Z$t  =  0.05  where  maximum  flame  temperature  occurs. 
Recently,  Seshadri  and  Peters  (1988)  performed  an  asymptotic  analysis  of  the  structure  of 
methane-air  diffusion  flames  using  realistic  flame  chemistry  to  show  that  the  peak 
temperature  occurs  at  a  value  of  Z  which  is  larger  than  Zst,  and  our  computations  support 
this  result  for  heptane-air  diffusion  flames. 

In  Figs.  8  and  9,  the  experimental  profiles  for  temperature  and  the  stable  species 
C7H16,  O2,  CO2,  H2O,  H2,  and  CO  are  compared  with  those  calculated  for  U  =  83.7  cm/s 
using  the  physical  variable  z  and  the  mixture  fraction  Z  as  the  independent  variable, 
respectively.  The  agreement  between  the  experimental  and  calculated  values  is  not  as  good, 
which  is  expected  because  of  greater  experimental  difficulties  in  spatial  resolution  at  this 
higher  strain  rate.  However,  the  agreement  is  qualitatively  similar  to  that  shown  in  Figs.  6 
and  7.  Fig.  10  shows  the  calculated  profiles  for  the  radicals  H,  OH,  O,  and  HO2  for  U  = 
30.5  cm/s  and  U  -  83.7  cm/s  and  is  discussed  below. 


3.2  Influence  of  Strain  Rate: 


To  evaluate  the  influence  of  strain  on  the  structure  of  heptane-air  flames,  some 
interesting  features  of  Figs.  5-10  are  shown  in  Table  2.  With  increasing  flow  velocity  of 
the  oxidizer  stream  (increasing  strain),  the  stagnation  plane  moves  toward  the  fuel  surface, 
and  the  flame  standoff  distance,  which  is  defined  as  the  distance  between  the  peak  value  of 
the  flame  temperature  and  the  fuel  surface,  decreases.  Also  with  increasing  strain,  the  fuel 
burning  rate  increases  because  the  higher  molar  flow  rate  >f  oxidizer  must  be  matched  with 
that  of  the  fuel,  such  that  at  the  flame,  the  flow  rate  of  the  reactants  are  in  stoichiometric 
proportions.  Table  2  shows  that,  with  increasing  strain,  the  maximum  flame  temperature 
decreases,  and  the  oxygen  concentration  at  that  position  increases.  Hence,  leakage  of 
oxidizer  through  the  reaction  zone  increases  with  strain.  However,  the  leakage  of  the  fuel, 
which  is  represented  by  the  fuel  concentration  at  maximum  temperature,  is  negligibly  small 
in  both  cases. 

These  results  are  consistent  with  those  of  asymptotic  analysis  done  by  Linan  (1974), 
Miller  et  al.  (1984),  Tangirala  et  al.  (1989),  Seshadri  and  Peters  (1988),  which  show  that 
with  increasing  strain,  the  flame  temperature  must  decrease,  and  the  leakage  of  the  reactants 
through  the  flame  must  increase.  However,  Linan  (1974)  showed  that,  when  the  chemistry 
for  the  hydrocarbon  flames  is  approximated  by  a  one-step,  irreversible  reaction,  the  leakage 
of  fuel  through  the  reaction  zone  is  considerably  higher  than  that  of  oxygen  with  increasing 
strain.  Whereas  results  using  realistic  chemistry,  as  reported  by  Seshadri  and  Peters 
(1988),  show  that  the  leakage  of  oxygen  through  the  reaction  zone  increases  with 
increasing  strain,  and  the  leakage  of  fuel  is  negligible.  Hence,  the  numerical  results  are 
qualitatively  consistent  with  asymptotic  analysis  results  using  realistic  chemistry. 
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Table  2  shows  that  the  value  of  Z  at  the  surface  of  the  liquid  pool  increases  with 
increasing  strain,  and  this  is  due  to  the  increasing  fuel  concentration  at  the  liquid  pool 
surface.  Table  2  also  oi.ows  considerable  discrepancy  between  the  measured  and  computed 
oxygen  concentration  at  maximum  temperature;  this  is  probably  due  to  inaccuracies 
associated  with  measuring  small  oxygen  concentrations  using  a  gas  chromatograph  since 
the  uncertainties  in  the  numerical  results  are  considerably  smaller.  This  conclusion  was 
also  reached  in  previous  comparisons  of  experimental  measurements  and  numerical 
calculations  for  methane-air  flames  done  by  Smooke  et  al.  (1986)  and  Puri  et  al.  (1987). 
Fig.  10  shows  that,  although  the  peak  concentration  of  the  radicals  H,  OH,  and  0  do  not 
change  much  with  increasing  strain,  the  peak  concentration  of  HO2  increases  sharply,  and 
this  is  consistent  with  the  steady-state  approximation  for  HO2,  as  reported  by  Seshadri  and 
Peters  (1988). 

3-3  Prospect  for  Asymptotic  Description: 

It  is  useful  to  compare  certain  aspects  of  the  structure  of  heptane-air  diffusion  flames 
relevant  for  asymptotic  description  with  similar  ones  of  methane-air  diffusion  flames.  If  as 
assumed  previously  for  methane-air  diffusion  flames,  the  concentration  of  all  species 
except  those  for  C7H16,  O2,  CO,  H2,  H,  H2O,  CO2,  and  N2  are  presumed  to  be  in  a 
dynamic  steady  state,  then  using  the  procedures  outlined  by  Peters  (1985),  it  can  be  shown 
that  the  chemical  kinetic  mechanism  for  heptane  oxidation  can  be  simplified  to  the  following 
four  overall  reactions: 


C7H16  +  10  H  +  7  H20  — >  7  CO  +  20  H2 

(1) 

CO  +  H20  <===>  CO2  +  h2 

(2) 

H  +  H  +  M  ->  H2  +  M 

(3) 

02  +  3  H2  <===>  2  H  +  2H20 

(4) 

%0 


The  rates  of  these  four  overall  reactions  can  be  related  to  the  rates  of  elementary  reactions 
shown  in  Table  1.  In  Fig.  8,  the  computed  profiles  for  the  species  appearing  in  the  above 
overall  reactions  are  shown  for  U  =  30.5  cm/s.  Two  separate  regions  are  identified  and 
labeled  as  the  "fuel  consumption  layer"  and  the  "oxygen  consumption  layer".  In  the  fuel- 
consumption  layer,  radicals  such  as  H  atoms  attack  the  fuel  to  form  CO  and  H2  in  accord 
with  the  overall  reaction  1.  Because  of  the  strong  affinity  of  the  fuel  to  the  radicals,  further 
oxidation  of  the  intermediate  species  H2  and  CO  is  inhibited  in  this  layer;  consequendy,  the 
concentrations  of  H2  and  CO  attain  a  maximum  value  in  the  fuel  consumption  layer.  After 
the  fuel  is  nearly  consumed  on  the  oxidizer  side  of  the  flame,  there  exists  an  "oxygen- 
consumption  layer",  where  the  overall  reactions  3-4  occur,  and  radicals  are  produced  to 
react  with  H2  and  CO  to  form  the  final  products  H2O  and  CO2.  Thus,  the  calculations 
clearly  show  that  certain  aspects  of  the  structure  of  heptane-air  diffusion  flame  relevant  for 
asymptotic  description  are  similar  to  those  of  methane-air  diffusion  flames,  as  shown  by 
Peters  (1985).  This  offers  considerable  simplification  because  it  implies  that  the  techniques 
developed  for  analyzing  the  structure  of  methane-air  flames  can  be  extended  to  heptane-air 
flames. 


4.  CONCLUSION 

Detailed  numerical  calculations  were  performed  to  determine  the  structure  of  heptane- 
air  diffusion  flames,  and  the  results  were  compared  with  experimental  measurements.  The 
following  remarks  summarize  the  major  conclusions: 

1 .  The  calculated  profiles  of  temperature  and  stable  species  were  found  to  agree  well 
with  the  experimental  measurements,  but  the  computed  peak  temperature  was  shifted 
slightly  and  was  attributed  to  the  inaccurate  description  of  the  flow  field.  However, 


better  alignment  was  obtained  when  the  profiles  were  plotted  using  the  mixture 
fraction  as  the  independent  variable. 

2.  With  increasing  strain,  the  burning  rate  increases,  the  peak  temperature  decreases, 
and  the  oxygen  leakage  increases. 

3 .  Key  aspects  of  the  structure  of  heptane-air  diffusion  flames  are  similar  to  those  of 
methane-air  diffusion  flames;  hence,  previously  developed  techniques  for 
asymptotic  analysis  can  be  extended  to  heptane-air  flames. 

The  research  summarized  above  has  been  submitted  for  publication  in  Combustion  Science 

and  Technology.  The  coauthor  of  the  publication  is  Ms.  Mary  Bui-Pham. 
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TABLE  1 


CHEMICAL  KINETIC  MECHANISM  FOR  HEPTANE  OXIDATION 


Reactions 

1.  H+02=0H+0 

2.  H2+0=0H+H 

3.  0H+H2=H20+H 

4.  0H+0H=H20+0 

5.  H2+M=H+H+M 

6.  H+0H+M=H20+M 

7.  H+02+M=H02+M 

8.  H+02+N2=H02+N2 

9.  H+H02=0H+0H 

10.  H+H02=H2+02 

11.  0+H02=0H+02 

12.  0H+H02=H20+02 

13.  H2+02=0H+0H 

14.  H+02+02=H02+02 

15.  C0+0H=C02+H 

16.  CH3+0=CH20+H 

17.  CH20+H=HC0+H2 

18.  CH20+M=HC0+H+M 

19.  CH20+0=HC0+0H 

20.  0H+CH20=HC0+H20 

21.  H+HCO=CO+H2 

22.  HC0+0=C0+0H 

23.  0H+HC0=C0+H20 

24.  02+HC0=C0+H02 

25.  HCO+M=CO+H+M 

26.  CH3+02=CH30+0 

27.  CH3+0H=CH20+H2 

28.  CH30+M=CH20+H+M 

29.  CH30+H=CH20+H2 

30.  CH30+0H=CH20+H20 

31.  CH30+0=CH20+0H 

32.  CH30+02=CH20+H02 

33.  CH3+02=CH20+0H 

34.  CH3HCO+H=CH3+CO+H2 

35.  CH3HC0+0=CH3+C0+0H 

36.  CH3HC0+0H=CH3+C0+H20 

37.  C3H6+0=CH3+CH3+C0 

38.  C3H6+OH=CH3HCO+CH3 

39.  C7H16+H=C7H15+H2 

40.  C7H16+0=C7H15+0H 

41.  C7H16+0H=C7H15+H20 

42.  C7H15=CH3+C3H6+C3H6 


A 

nj 

E 

2.20e+14 

0. 

16800 

1.80e+10 

1.00 

8826 

1.17e+09 

1.30 

3626 

6.00e+08 

1.30 

0 

2.23e+12 

0.50 

92600 

7.50e+23 

-2.60 

0 

2.10e+18 

-1.00 

0 

6.70e+19 

-1.42 

0 

2.50e+14 

0. 

1900 

2.50e+13 

0. 

700 

4.80e+13 

0. 

1000 

5.00e+13 

0. 

1000 

1.70e+13 

0. 

47780 

6.70e+19 

-1.42 

0 

1.51e+07 

1.30 

-758 

6.80e+13 

0. 

0 

3.31e+14 

0. 

10500 

3.31e+16 

0. 

81000 

1.81e+13 

0. 

3082 

7.53e+12 

0. 

167 

4.00e+13 

0. 

0 

1.00e+12 

0. 

0 

5.00e+12 

0. 

0 

3.00e+12 

0. 

0 

1.60e+14 

0. 

14700 

7.00e+12 

0. 

25652 

7.50e+12 

0. 

0 

2.40e+13 

0. 

28812 

2.00e+13 

0. 

0 

1.00e+13 

0. 

0 

1.00e+13 

0. 

0 

6.30e+10 

0. 

2600 

5.20e+13 

0. 

34574 

4.00e+13 

0. 

4200 

5.00e+12 

0. 

1900 

1.00e+13 

0. 

0 

5.00e+12 

0. 

454 

1.00e+13 

0. 

0 

6.10e+14 

0. 

8469 

1.60e+14 

0. 

4569 

1.70e+13 

0. 

957 

3.70e+13 

0. 

28708 

Tt- 


43. 

02+M=0+0+M 

70. 

C2H2+OH=CH2CO+H 

44. 

CH3+H=CH4 

71. 

C2H2+0=CH2+C0 

45. 

CH4+H=CH3+H2 

72. 

CH2CO+H=CH3+CO 

46. 

CH4+0=CH3+0H 

73. 

CH2C0+0=HC0+HC0 

47. 

CH4+0H=CH3+H20 

74. 

CH2C0+0H=CH20+HC0 

48. 

CH4+02=CH3+H02 

75. 

CH2CO+M=CH2+CO+M 

49. 

C0+0+M=C02+M 

76. 

C2H2+0=C2H0+H 

50. 

C0+02=C02+0 

77. 

C2H0+0=C0+C0+H 

51. 

CH2+H=CH+H2 

78. 

C2HO+H=CH2+CO 

52. 

CH2+0=C0+H+H 

79. 

C2H2+0H=C2H+H20 

53. 

CH2+02=C02+H+H 

80. 

C2H+0=C0+CH 

54. 

CH2+CH3=C2H4+H 

81. 

C2H+H2=C2H2+H 

55. 

CH+0=C0+H 

82. 

C2H+02=C0+HC0 

56. 

CH+02=C0+0H 

83. 

C2H5+CH3=C3H8 

57. 

CH3+CH3=C2H6 

84. 

C3H8+H=N*C3H7+H2 

58. 

CH3+CH3=C2H5+H 

85. 

C3H8+H=I*C3H7+H2 

59. 

C2H6+H=C2H5+H2 

86. 

C3H8+0=N*'C3H7+0H 

60. 

C2H6+0=C2H5+OH 

87. 

C3H8+0=I*C3H7+0H 

61. 

C2H6+0H=C2H5+H20 

88. 

C3H8+0H=N*C3H7+H20 

62. 

C2H5+02=C2H4+H02 

89. 

C3H8+0H=I*C3H7+H20 

63. 

C2H5=C2H4+H 

90. 

N*C3H7+H=C3H8 

64. 

C2H4+0=HC0+CH3 

91. 

I*C3H7+H=C3H8 

65. 

C2H4+0H=C2H3+H20 

92. 

N*C3H7+02=C3H6+H02 

66. 

C2H4+H=C2H3+H2 

93. 

I*C3H7+02=C3H6fH02 

67. 

C2H3+H=C2H2+H2 

94. 

N*C3H7=C2H4+CH3 

68. 

C2H3+02=C2H2+H02 

95. 

N*C3H7=C3H6+H 

69. 

C2H3=C2H2+H 

96. 

I*C3H7=C3H6+H 

Rate  coefficients  are  in  the  form  kj  =  Aj  T"*  exp  (-  Ej  /  RT).  Units  are  moles, 
cubic  centimeters,  seconds,  Kelvins  and  calories/mole. 
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Fig.  3  Comparison  between  concentration  profiles  for  H,  O,  OH,  and  HO2  calculated  using  reactions  1- 
96  (dashed  lines)  with  those  calculated  using  reactions  1-42  (solid  lines)  for  a=9.84  s-* 
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Fig.  4  Schematic  illustration  of  the  principal  path  of  oxidation  of  heptane 


Fig.  5  Calculated  profile  of  the  axial  component  of  the  velocity  for  YC>2«  =  0.185,  U=30.5  cm/s, 
U=83.7  cm./s 
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Fig.  6.  Comparison  between  experimental  measurements  (shown  as  points)  and  computed  values  (lines) 
for  profiles  of  temperature  and  stable  species:  C7H16,  O2,  CO2,  H2,  CO,  and  H2O  for 
YO->  -fURS  l  i-ro  /o 
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Fig.  7.  Comparison  between  experimental  measurements  (shown  as  points)  and  computed  values  (line 
for  profiles  of  temperature  and  stable  species:  C7Hj6,  O2,  CO2,  H2,  CO,  and  H2O  using  the 
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Fig.  8.  Comparison  between  experimental  measurements  (shown  as  points)  and  computed  values  (lines) 
for  profiles  of  temperature  and  stable  species:  C7H16,  O2,  CO2,  H2,  CO,  and  H2O  for 
YG2«=0.185,  U=83.7  cm/s 


Fig.  9.  Comparison  between  experimental  measurements  (shown  as  points)  and  computed  values  (lines) 
for  profiles  of  temperature  and  stable  species:  C7H16,  O2,  CO2,  H2,  CO,  and  H2O  using  the 
mixture  fraction  Z  as  the  independent  variable,  for  YO2oo=0.185,  U=83.7  cm/s 


Fig.  10.  Calculated  profiles  for  the  radicals  H,  O,  OH,  and  HO2  forY02~  =  0.185,  U=30.5  cm/s, 
U=83.7  cm./s 
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Table  2  Certain  aspects  illustrating  the  influence  of  strain  on  structure 
of  heptane-air  diffusion  flames 


CHAPTER  III 


ANALYSIS  OF  THE  STRUCTURE  AND  MECHANISMS  OF  EXTINCTION  OF 
COUNTERFLOW  METHANOL-/.  DIFFUSION  FLAME 

Numerical  calculations  were  performed  to  determine  the  structure  and  to 
clarify  the  mechanisms  of  extinction  of  diffusion  flames  stabilized  between 
counterflowing  streams  of  methanol  and  air.  The  calculations  were  performed 
at  a  value  of  the  thermodynamic  pressure  equal  to  one  atmosphere,  with 
different  values  of  the  rate  of  strain  ,  and  with  two  different  chemical  kinetic 
mechanisms.  We  will  refer  to  these  two  mechanisms  as  "mechanism  a"  and 
"mechanism  b".  Mechanism  a  and  mechanism  b  have  the  same  set  of 
elementary  reactions,  but  the  rate  constants  for  these  elementary  reactions 
were  obtained  from  two  different  references.  Temperature  profiles, 
concentration  profiles  of  various  species,  rates  of  production  and  destruction 

of  various  species,  and  rates  of  various  reactions  were  plotted  as  a  function  of 
the  axial  coordinate  using  mechanism  a  and  mechanism  b.  Both  chemical 
mechanisms  show  that  the  structure  can  be  subdivided  into  three  regions,  the 
fuel  consumption  region  where  the  reaction  proceeds  via  the  path  CH3OH  - 

CH2OH  -  CH2O  -  CO,H2;  the  H2  -  CO  oxidation  region  where  the  compounds  H2 
and  CO  oxidize  to  form  H2O  and  CO2;  and  the  radical  destruction  region  where 
radicals  are  destroyed  by  three  body  reactions. 

If  mechanism  a  is  used,  we  conclude  that  at  low  rates  of  strain  the 

concentration  of  CH2OH  and  HCO  arc  in  steady-state  and  in  partial  equilibrium 

is  assumed  for  certain  reactions  there  exists  algebraic  relations  among  the 
concentrations  of  the  radicals  OH,  H,  and  0.  As  the  rate  of  strain  is  increased, 
HCO  is  no  longer  in  steady-state  and  no  solution  was  obtained  at  a  value  of  the 
strain  rate  greater  than  521  s*  *.  However,  if  mechanism  b  is  used,  the 
concentration  of  only  HCO  is  in  steady-stale,  and  there  also  exist  algebraic 
relations  among  the  concentrations  of  the  radicals  OH,  H,  and  0.  As  the  rate  of 
strain  is  increased,  no  solution  was  obtained  for  a  value  of  the  strain  rale 
greater  than  168  s' ^ ,  and  we  speculate  that  extinction  of  the  flame  is  caused 
due  to  a  large  value  of  the  activation  energy  for  a  reaction  controlling  the 
pyrolysis  of  CH2OH  to  CItoO. 


The  research  summarized  above  has  been  published  in  Combustion  and 
Flame  76,  pp  111-132,  1989.  The  coauthors  of  the  publication  are  Professor  C. 
Trevino,  and  Professor  M.  D.  Smooke 


CHAPTER  IV 


ASYMPTOTIC  STRUCTURE  AND  EXTINCTION  OF  METHANE-AIR 

DIFFUSION  FLAMES 

The  asymptotic  structure  of  a  counterflow  methane-air  flame  is  analyzed 
using  a  three-step  chemical  kinetic  mechanism, 

I  CH4  +  02  -->  CO  +  H2  +  H2O 

II  CO  +  H2O  <=>  CO2  +  H2 

III  O2  +  H2  «>2H20 

which  was  deduced  in  a  systematic  way  through  steady  state  and  partial 
equilibrium  assumptions  from  a  detailed  chemical  kinetic  mechanism  for 
oxidation  of  methane.  The  rates  for  the  three  steps  are  related  to  the  rates  of 
the  elementary  reactions  .  The  outer  structure  of  the  diffusion  flame  is  the 
classical  Burkc-Schumann  structure  governed  by  the  overall  one-step 
reaction  CH4  +  202  -■>  CO2  +  2H2O,  with  the  flame  sheet  positioned  at  Z  =  Zst, 
where  Z  is  the  mixture  fraction  used  as  the  independent  variable  in  the 
analysis.  The  inner  structure  consists  of  a  thin  H2  -  CO  oxidation  layer  of 
thickness  0(E)  toward  the  lean  side,  a  thin  nonequilibrium  layer  for  the  water 
gas  shift  reaction  of  thickness  of  0(v),  and  a  thin  fuel  consumption  layer  of 
thickness  0(8)  toward  the  rich  side.  These  layers  result,  respectively,  in  the 
limit  of  large  values  for  the  Damkohlcr  number  characterizing  the  rate  of 
reaction  III,  III,  and  I,  while  the  ratios  of  the  activation  temperature  to  the  gas 
temperature  for  the  three  reactions  arc  assumed  to  be  of  order  unity.  We  also 
find  that  E  >  V  >  8.  The  results  of  the  asymptotic  analysis  yield  values  of  the 

temperature  and  the  main  species  at  the  fuel  consumption  layer  as  a  function 
of  the  scalar  dissipation  rate  %st.  wc  therefore  obtain  the  upper  branch  f  the 

classical  S-shapcd  curve  where  the  maximum  flame  temperature  is  plotted  as  a 
function  of  %S['L  The  scalar  dissipation  rate  at  quenching  %q  is  derived  from 

the  S-shapcd  plot  and  its  relation  to  the  laminar  burning  velocity  is  discussed. 
A  comparison  of  the  diffusion  flame  structure  with  that  of  a  premixed  flame 
shows  that  the  rich  part  of  the  diffusion  flame  corresponds  to  the  upstream 
part  of  the  premixed  flame  while  its  lean  part  corresponds  to  the  downstream 
part.  First  the  kinetic  scheme  is  based  on  the  most  important  (principal) 
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reactions  to  derive  the  basic  structure  .  When  a  number  of  additional 
elementary  chemical  reactions  are  added  the  results  of  the  asymptotic  analysis 
are  found  to  be  in  very  good  agreement  with  previous  numerical  calculations 
that  used  a  complete  kinetic  mechanism,  as  well  as  with  experiments. 

The  research  summarized  above  has  been  published  in  Combustion  and 
Flame  73,  pp  23  -  44,  1988.  The  coauthor  of  this  publication  is  Professor  N. 
Peters.  The  international  collaboration  was  supported  by  the  National  Science 
Foundation  Grant  number  NSF-INT-86-09939 


CHAPTER  V.  METHANE  -  AIR  PREMIXED  FLAMES 


V.  1  COMPARISON  BETWEEN  EXPERIMENTAL  MEASUREMENTS  AND 
NUMERICAL  CALCULATIONS  OF  THE  STRUCTURE  OF  COUNTERFLOW, 
DILUTED,  METHANE-AIR,  PREMIXED  FLAMES. 

ABSTRACT 

In  this  chapter  results  of  experimental  and  numerical  study  of  the 
structure  of  premixed  flame  is  reported.  The  fuel  used  in  this  study  is  gaseous 
methane.  Experiments  were  performed  in  premixed  flames  stabilized  between 
counterflowing  streams  of  a  combustible  mixture  of  methane  and  air  diluted 
with  nitrogen  and  a  stream  which  consists  of  reaction  products  from  another 
premixed  methane-air  flame.  Critical  conditions  of  extinction  of  this  flame 
were  measured  over  a  wide  parametric  range.  The  detailed  structure  of  the 
flame  was  characterized  by  measuring  the  concentration  profiles  of  various 
stable  species  by  removing  gas  samples  from  the  flame  using  a  quartz 
microprobe.  The  gas  samples  were  then  analyzed  employing  gas 
chromatographic  techniques.  Temperature  profiles  were  measured  by  use  of 
coated  thermocouples.  The  measurements  were  corrected  by  radiative  losses 
from  the  thermocouple  beads.  Detailed  numerical  calculations  were  performed 
to  determine  the  structure  of  the  flame  and  the  results  were  compared  with 
experimental  methods. 

It  is  shown  that  the  basic  structure  of  most  hydrocarbon  flames  including 
that  of  methane  arc  similar.  Hence,  some  of  the  conclusions  deduced  here  can 
be  extended  to  other  premixed  flames  burning  more  complex  hydrocarbon 
fuels.  Methane  was  chosen  for  the  present  study  because  gaseous  fuels  permit 
experimental  measurements  to  be  made  over  a  wide  parametric  range. 


A  manuscript  describing  the  research  outlined  in  this  chapter  will 
appear  in  the  proceedings  of  the  Twenty-Third  (International)  Symposium  on 
Combustion,  The  Combustion  Institute  1990.  The  coauthors  of  this  manuscript 
are  Professor  M.  D.  Smooke,  Captain  J.  Crump,  and  Dr.  V.  Giovangigli.  The 
authors  acknowledge  the  invaluable  assistance  of  Dr.  Venkat  Tangirala  with 
the  gas  chromatographic  measurements.  The  authors  also  acknowledge  the 
helpful  discussions  which  they  had  with  Professors  F.  A.  Williams,  N.  Peters, 
and  P.  A.  Libby.  The  research  performed  at  Yale  University  was  supported  by 
the  United  States  Department  of  Energy  and  the  Air  Force  Office  of  Scientific 
Research. 


V.2  THE  INNER  STRUCTURE  OF  METHANE-AIR  FLAMES 


The  inner  structure  of  a  methane-air  premixed  flame  is  analyzed  using  a  re 


duccd  four-step  chemical  kinetic  mechanism 

CH4  +  2H  +  H20  — >  CO  +  4H2  I 

CO  +  h2o  <— >  co2  +  h2  II 

H  +  H  +  M  — >  H2  +  M  III 

02  +  3H2  <==>  2H20  +  2H  IV 


The  rates  of  these  four  steps  are  related  to  the  rates  of  elementary  reactions 
appearing  in  the  Cj  chain  mechanism  for  oxidation  of  methane.  The  inner  layer  is 

thin  with  reactions  I  - 1 V  occurring  in  this  layer,  and  is  embedded  between  a  chemi¬ 
cally  inert  upstream  layer  and  a  broader  (  but  asymptotically  thin  )  downstream 
layer  where  reactions  II,  III  and  IV  occur  and  H2  and  CO  are  oxidized. 


The  analysis  reported  here  extends  a  previous  analysis  by  Peters  and  Williams 
of  the  structure  of  premixed  methane  -  air  flames,  where  a  reduced  three  -  step 
chemical  kinetic  mechanism  was  used.  In  the  equations  describing  the  structure  of 

the  inner  layer  a  parameter  0)  appears  which  represents  the  ratio  of  the  thickness  of 
the  fuel  consumption  layer  to  the  thickness  of  the  radical  consumption  layer  of  the 
previous  analysis  by  Peters  and  Williams.  Analytical  solutions  for  the  burning  veloc¬ 
ity  eigenvalue  L  are  obtained  in  the  limit  0)  -->  0  and  CO  -->  infinity,  and  by  use  of  nu¬ 
merical  integration  an  approximation  for  L  is  obtained  as  a  function  of  CO,  which  in¬ 


cludes  limiting  expressions  for  CO  -->  0  and  CO  -->  infinity.  The  expressions  for  L  con¬ 
tains  a  number  of  parameters,  which  represent  the  influence  of  a  number  of  ele¬ 


mentary  chemical  reactions.  In  particular,  a  parameter  defined  as  (i  in  the  analysis  is 
found  to  have  a  significant  influence  on  the  value  of  L  and  consequently  on  the 
burning  velocity,  and  the  influence  of  this  quantity  increases  with  increasing  pres¬ 


sure.  The  parameter  |1  represents  the  influence  of  the  backward  steps  of  the  reac¬ 
tions  CH4  +H<==>CH3  +  H2  and  CH4+OH<=>CH3  +  H-,0. 


Using  the  results  of  the  analysis  the  burning  velocity  was  calculated  for  a  stoi¬ 
chiometric  methane-  air  flame  for  values  of  the  pressure  p  between  1  atm  and  80  atm. 
at  p  =  1  atm  the  calculated  burning  velocity  was  38  cm/s  in  good  agreement  with  ex¬ 
perimental  results.  The  burning  velocity  decreased  with  increasing  pressure  again 
in  agreement  with  experimental  measurements. 

The  research  summarized  above  has  been  published  in  Combustion  and  Flame 
81,  pp  96-118,  1990.  The  coauthor  of  this  publication  is  Professor  Norbert  Peters. 


V.3  THE  ASYMPTOTIC  STRUCTURE  OF 
NONSTOICHIOMETRIC  METHANE-AIR  FLAMES 


Abstract 

The  asymptotic  structure  of  methane-air  premixed  flames  is  analyzed  for  all  values  of 
the  equivalence  ratio,  using  a  reduced  chemical-kinetic  mechanism  that  can  be  represented 
as 

I  CH4  +  2H  +  H2O  CO  +  4H2 

II  CO  +  H20  2  C02  +  H2 

m  h  +  h  +  m->h2  +  m 

IV  02  +  3H2  £  2H20  +  2H 

The  rates  of  these  four  steps  are  related  to  the  rates  of  elementary  reactions  appearing  in  the 
Ci-chain  mechanism  for  oxidation  of  methane.  In  the  analysis  the  overall  flame  structure  is 
subdivided  into  four  zones,  a  preheat  zone,  an  inner  layer,  a  H2-CO  oxidation  layer  and  a 
post-flame  zone.  The  inner  layer  is  located  between  the  preheat  zone  and  the  oxidation 
layer.  The  inner  layer  is  thin  in  comparison  to  the  preheat  zone,  and  in  this  layer  finite-rate 
reactions  related  to  the  consumption  of  the  fuel  are  considered.  Downstream  from  this 
inner  layer  there  exists  a  broader  (but  asymptotically  thin)  oxidation  layer  where  reaction  II, 
III  and  IV  occur  and  H2  and  CO  are  oxidized  to  form  H20  and  CO2. 

The  asymptotic  analysis  employed  for  the  inner  layer  is  essentially  identical  to  the 
previous  analysis  of  Seshadri  and  Peters.  However,  the  analysis  of  the  H2-CO  oxidation 
layer  is  different  because  the  values  of  the  concentrations  of  the  species  H2,  CO  and  02  are 
not  necessarily  negligibly  small  in  the  post-flame  zone  downstream  from  the  oxidation 
layer.  The  results  of  the  asymptotic  analysis  are  used  to  calculate  the  burning  velocity  of 
the  flame  as  a  function  of  the  equivalence  ratio  for  various  values  of  the  pressure.  An 
apparent  overall  activation  energy  is  deduced  from  the  analysis,  and  the  corresponding 
Zel'dovich  number  increases  rapidly  as  the  equivalence  ratio  departs  from  unity.  The  large 
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value  of  the  Zel'dovich  number  for  very  fuel-rich  and  very  fuel-lean  flames  makes  these 
flames  much  more  sensitive  to  heat  loss  or  flame  stretch.  This  result,  coupled  with  the 
finding  from  the  analysis  that  a  fundamental  chemical-kinetic  mechanism  for  extinction 
does  not  exist  at  the  experimentally  observed  flammability  limits,  suggests  that  flammability 
limits  are  due  to  heat  losses  for  methane-air  flames. 

1 .  Introduction 

The  asymptotic  structure  of  stoichiometric,  laminar,  premixed  methane-air  flames  has 
been  investigated  previously  [1,2]  using  a  reduced  chemical-kinetic  mechanism  to  describe 
the  gas-phase  oxidation  of  this  fuel.  This  reduced  mechanism  was  deduced  systematically 
from  a  detailed  chemical-kinetic  mechanism  by  assuming  that  the  concentrations  of  a 
number  of  intermediate  species  are  in  a  dynamic  steady  state  [3].  In  the  analysis  of  Peters 
and  Williams  [1],  a  three-step  mechanism  was  employed,  and  an  asymptotic  analysis  was 
performed  using  only  the.  principal  elementary  reactions  to  characterize  the  overall  rates  of 
the  reduced  mechanism.  Results  for  the  burning  velocity  as  a  function  of  pressure  and 
orders  of  magnitude  of  the  concentrations  of  all  intermediates  were  obtained  [1].  The 
results  of  the  asymptotic  analysis  were  improved  by  including  a  number  of  additional 
chemical  reactions  in  the  overall  rates  of  the  reduced  mechanism  through  numerical 
evaluation  of  an  integral  involving  an  iteration  for  determination  of  the  integrand.  In  the 
analysis  of  Seshadri  and  Peters  [2],  a  four-step  mechanism  was  used,  and  the  effects  of 
additional  reactions  were  represented  by  additional  algebraic  parameters  to  avoid  the 
iteration  of  Ref.  1,  which  would  be  inapplicable  with  the  four-step  mechanism.  Results 
were  obtained  for  the  burning  velocity  of  the  flame  for  values  of  pressure  p  including  the 
range  1  atm  to  40  atm,  addressed  here.  At  p  =  1  atm  the  calculated  burning  velocity  was  35 
cm/s,  in  reasonable  agreement  with  experimental  results  [2].  The  burning  velocity 
decreased  with  increasing  pressure,  again  in  agreement  with  experimental  measurements. 

Here  the  analyses  reported  in  Refs.  1  and  2  are  extended  to  all  values  of  the 
equivalence  ratio  <|).  Only  terms  of  leading  order  in  the  asymptotics  are  retained  because  the 
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contributions  of  the  higher-order  terms  have  been  estimated  to  be  20%  or  less,  comparable 
with  inaccuracies  introduced  by  other  approximations  in  the  analysis.  The  intent  is  to 
illustrate  the  character  of  the  analysis  and  of  the  results  and  to  provide  an  indication  of  the 
errors  that  may  be  expected  from  the  asymptotics  at  leading  order.  Investigation  of  the 
greater  accuracy  that  can  be  achieved  by  proceeding  to  higher  order  is  not  an  objective  of 
the  present  paper.  Owing  to  the  neglect  of  C2  chemistry,  the  results  of  the  analysis  are 
expected  to  be  notably  inaccurate  for  fuel-rich  flames,  and  therefore  potential  improvements 
in  the  asymptotics  for  fuel-rich  conditions  are  not  thoroughly  explored.  The  asymptotic 
analysis  and  notation  used  here  closely  follow  the  development  in  Ref.  2. 

2.  Structure  of  the  Post-Flame  Zone 

Downstream  of  the  flame,  in  the  post-flame  zone  the  reaction  products  are  in  chemical 
equilibrium,  and  the  temperature  is  equal  to  the  adiabatic  flame  temperature.  The 
equilibrium  concentrations  of  the  products  are  determined  as  described  elsewhere  [4]  by 
assuming  that  only  the  species  O2,  H2,  CO,  CO2  and  H2O  are  present  and  are  in  chemical 
equilibrium,  and  that  the  enthalpy  and  the  element  mass  fractions  are  equal  to  those  in  the 
unbumt  gas.  For  given  values  of  <j),  p  and  the  initial  temperature,  Tu,  the  adiabatic  flame 
temperature,  Tb,  and  the  mass  fraction  of  species  i  in  the  post-flame  zone,  Yjb,  can  be 
calculated.  The  equivalence  ratio  is  related  to  the  mixture  fraction  Z  by 


(1-Zst) 
(1-Z)  ’ 


0) 


where  Zst  =  0.055  is  the  stoichiometric  mixture  fraction.  In  Fig.  1  results  of 
thermochemical  calculations  are  plotted  showing  the  conditions  in  the  post-flame  zone  as 
functions  of  <!>. 

Shown  in  Fig.  la,  in  addition  to  Tb,  is  the  adiabatic  flame  temperature  Tc  for 
complete  combustion  with  Yh2  =  Yco  =  0,  which  corresponds  to  the  overall  reaction 
CH4  +  2O2  — >  CO2  +  2H2O.  In  this  limit  the  products  would  contain  unbumt  O2  for  fuel- 
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lean  flames  and  unbumt  CH4  for  fuel-rich  flames.  If  the  heat  of  combustion  for  this  overall 
reaction  per  mole  of  CH4  consumed  is  denoted  by  (-AH),  then 

dT  =  |  (-AH)Yfu/Wf  for  4>  <  1 

P,P  |  (-AH)Yo2U/(2W0)  for  <j>  >  1 

where  cP(p  is  the  specific  heat  at  constant  pressure  of  the  product  mixture  (N2,  CO2,  H2O 
and  O2  or  CH4),  Ypu  and  Yq^  denote  the  mass  fractions  of  fuel  and  oxidizer,  respectively, 

in  the  initial  reactant  stream,  and  Wj  is  the  molecular  weight  of  species  i.  Figure  1  shows 
that  the  discontinuity  in  slope  at  <j)  =  1,  predicted  by  this  formula,  is  smoothed  when  the 
equilibrium  chemistry  is  considered,  and  there  is  a  region  in  which  the  products  contain  O2, 
H2  and  CO  simultaneously.  Expansions  later  will  be  performed  about  Tc  rather  than  Tb. 

3.  Reduced  Chemical-Kinetic  Mechanism 

A  chemical-kinetic  mechanism  describing  the  oxidation  of  methane  is  shown  in 
Table  1.  The  rate  parameters  here  differ  somewhat  from  those  employed  in  Ref.  1  but  are 
the  same  as  those  of  Ref.  2  and  are  obtained  from  a  recent  [5]  reevaluation  of  the  rate  data. 
Employing  the  procedure  described  in  Ref.  6,  a  four-step  mechanism  can  be  deduced  from 
this  mechanism  by  assuming  that  there  exists  a  dynamic  steady  state  for  the  concentrations 
of  the  species  O,  OH,  HO2,  H2O2,  CH3,  CH2O,  CH3O  and  HCO.  The  four-step 
mechanism  can  be  written  as 

I  CH4  +  2H  +  H2O  £  CO  +  4H2, 

II  CO  +  H2O  £  CO2  +  H2, 

III  H  +  H  +  M-»H2  +  M, 

IV  02  +  3H2  £  2H20  +  2H. 
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The  reaction  rates  Wk,  k  =  I,  n,  IE,  IV  for  the  overall  reactions  I  -  IV  can  be  related 
to  the  reaction  rates  of  the  elementary  reactions  wn,  n  =  1,2, .. . ,  25  shown  in  Table  1  and 
are 

Wl  =  Wi0f-Wi0b  +  Wiif-Wiib  +  Wi2f-Wi2b, 
wn  =  w9f-w9b, 

Wm  =  W5  -  WlOf  +  WlOb  +  W|6  -  W18  +  W19  -  W22f  +  W22b  +  W24  +  W25 , 

Wiv  =  Wif-Wib  +  W6  +  Wi8  +  W22f-W22b, 


where  the  subscripts  f  and  b  identify  forward  and  backward  rates,  respectively.  The  values 
of  the  reaction  rates  wn  shown  in  Eq.  (2)  are  proportional  to  the  product  of  the 
concentration  of  the  reactants  and  the  rate  constant  kn  of  the  elementary  reaction.  Results  of 
numerical  calculations  have  shown  that  reactions  lOf,  19, 21, 22f,  22b,  23f  and  23b  have 
only  a  minor  influence  on  the  burning  velocity  of  premixed  flames;  therefore,  as  in  Ref.  2, 
they  are  neglected  in  the  analysis  reported  here.  In  addition  to  the  steady-state 
assumptions,  as  in  previous  analyses  [1,2],  the  elementary  reaction  3  shown  in  Table  1  is 
assumed  to  be  in  equilibrium  yielding  the  algebraic  relation  Coh  =  YCh,  where  Q  is  the 
molar  concentration  of  species  i  and 


3  CH2O 
K3Ch2 


(3) 


in  which  K3  =  0.23  exp  (7530/T)  is  the  equilibrium  constant  of  the  elementary  reaction  3. 
4.  The  Conservation  Equations  for  a  Steady  Premixed  Flame 

For  a  steady,  planar,  adiabatic  deflagration  at  low  Mach  number  the  equation  of 
motion  implies  that  the  pressure  is  essentially  constant.  The  equation  for  mass 
conservation  can  be  written  as 

pv  =  puvu,  (4) 

where  p  is  the  density  and  v  the  gas  velocity.  The  index  u  denotes  conditions  in  the 
unbumt  gas.  Lewis  number.;  for  species  i  are  defined  as  L;  =  A./(pcpDj),  where  X  is  the 
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thermal  conductivity  and  cp  is  the  mean  specific  heat;  the  diffusion  coefficient  D,  is  taken  to 
be  that  of  species  i  with  respect  to  nitrogen,  and  the  binary-diffusion  approximation  is 
employed.  The  values  of  the  Lewis  numbers  for  all  species  are  assumed  to  be  constant. 
Using  the  notation  of  the  previous  analysis  [1,2],  the  non-dimensionalized  species  and 
energy  balance  equations  can  be  written  as 

Lp(Xp)  =  -  coi, 

Lh(Xh)  =  -  2o)i  -  2cuni  +  2o)rv , 

^h2(^Hi)  =  4q>i  +  con  +  com  -  3o)iv » 

=  -  coi  -  on  +  2corv , 

L<h.(  XqJ  =  -corv, 

^co(Xco)  =  toi-con, 
icOjfXcOi)  =  con, 

U(x)  =  Qi(0i  +  Qn(0ii  +  Qffltiim +Qiy®vi* 


The  operators  are  defined  as  Lj  s  d/dx  -  (l/Lj)d2/dx2,  and  L  =  d/dx  -  d2/dx2.  The  non- 
dimensional  independent  variable  x  is  related  to  the  spatial  coordinate  x'  as 
x  =  I  (pvcp/A,)dx',  and  the  quantities  Xj  and  x  are  related  to  the  mass  fraction  of  species 

Jo 

i,  Yi  and  the  gas  temperature,  T  as 


x 


T-Tu 

Tc-V 


(6) 


where  the  subscript  F  denotes  the  fuel.  In  the  analysis  the  average  molecular  weight  W  is 
assumed  to  be  a  constant  equal  to  27.62  kg/kmol;  hence  Xi  is  the  mole  fraction  of  species  i 
divided  by  the  initial  mole  fraction  of  the  fuel.  The  non-dimensionalized  reaction  rates  C0k 


and  the  non-dimensionalized  heats  of  reaction  Qk  of  the  reduced  four-step  mechanism  are 
defined  as 


_  kWFwk  q  _  YFU(-AHk) 

CpYFu(puVu)2  ’  k  cp(Tc  -  Tu)  Wp ' 


7 


The  non-dimensionalization  of  Eq.  (7)  will  also  be  applied  to  the  rates  of  the  elementary 
steps.  Since  assuming  steady  states  and  negligible  concentrations  for  CO,  H2  and  H 
enables  the  overall  reaction  CH4  +  2O2  — >  CO2  +  2H2O  to  be  deduced  from  the  four-step 
mechanism  by  adding  twice  reaction  IV  to  the  sum  of  reactions  I,  II  and  HI,  it  follows  from 
the  definition  of  Qk  given  in  Eq.  (7)  that  Qi  +  Qn  +  Qin  +  2Qrv  =  1. 

A  schematic  illustration  of  the  presumed  structure  of  the  premixed  flame  is  shown 
in  Fig.  2.  It  consists  of  a  chemically  inert  preheat  zone,  followed  by  the  inner  layer,  the 
H2-CO  oxidation  layer  and  the  equilibrium  post-flame  zone.  The  structures  of  the  inner 
layer  and  the  oxidation  layer  are  similar  to  those  shown  in  Ref.  2.  In  the  inner  layer  all  the 
hydrocarbon  chemistry  occurs,  resulting  in  the  formation  of  H2  and  CO  as  well  as  some 
H2O  and  CO2.  In  the  oxidation  layer  H2  and  CO  are  oxidized  to  form  H2O  and  CO2.  The 
concentration  of  fuel  is  zero  in  the  oxidation  layer,  and  the  H-radicals  are  in  steady  state.  In 
the  post-flame  zone  downstream  of  the  H2-CO  oxidation  layer,  H2  and  CO  are  in  partial 
equilibrium  according  to  reaction  n.  However,  there  exists  a  sublayer  between  the  inner 
layer  and  the  oxidation  layer  where  reaction  II  is  not  in  equilibrium. 

5 .  Asymptotic  Analysis  of  the  Inner  Layer 

Since  the  asymptotic  analysis  of  the  inner  layer  is  identical  to  that  of  Ref.  2,  only  the 
results  of  the  previous  analysis  will  be  shown  here.  The  thickness  of  this  layer  is 
presumed  to  be  of  order  8,  where 


5  = 


kifki3Xp2 

t  / 

knfki3Lp 


(8) 


in  which  the  superscript  0  implies  that  these  quantities  are  evaluated  at  the  origin,  x  =  0, 
which  is  taken  to  coincide  with  the  location  of  the  inner  layer,  and  at  this  point  T  =  T°  and 
X;  =  xf.  Other  quantities  in  Eq.  8  are 
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k13  ~  k13  +  i(k2fk10b  +  ^1  lb)/^lf]  (^Ha/  Od2)  » 

kllf  H  kllf  +  ,Yki2f’ 
k13  =  k13+Yki8/Ki, 
k2f  =  k2f  +  YK3k4b  , 
kllb  =  kllb +  ^K3k12b' 

The  value  of  Xq2  is  of  order  e,  where  e  is  a  measure  of  the  thickness  of  the  H2-CO 

»  „  , 

oxidation  layer.  In  addition,  the  ratio  klf/kllf  is  small,  and  k13/k13  is  roughly  of  order 
unity;  hence  8  will  be  presumed  to  be  smaller  than  e,  and  the  ordering  8  «  e  used  in  the 
previous  analyses  [1,2]  is  retained  here. 

Following  the  analysis  in  Ref.  2  a  quantity  L  which  contains  the  burning  velocity  is 
defined  as 

L  =  A82(kjlfk13/kj3)°LFR0  ,  (10) 

where 

R  a 

The  concentration  of  the  third  body  Cm  appearing  in  Eq.  (1 1)  can  be  written  in  terms  of  the 

N 

chaperon  efficiency  iqi  of  species  i  (i  =  1, . . .  N)  as  Cm  =  [  pW/(RTj]  ^  where 

i  =  l 

✓N 

the  gas  constant  is  R  =  82.05  atmcm3/(mol  K),  and  p  is  in  atmospheres.  In  the  previous 

analysis  [2]  two  limiting  structures  for  the  inner  layer  were  identified  depending  on  the 
value  of  a  parameter  co  defined  as 

CO  =  28LH/R°.  (12) 
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The  quantity  (0  represents  the  ratio  of  the  thickness  of  the  fuel-consumption  layer  (of  order 
8),  where  the  reaction  I  is  presumed  to  occur,  to  the  thickness  of  the  radical-consumption 
layer  of  the  previous  analysis  [1]. 

A  schematic  illustration  of  the  presumed  structure  of  the  inner  layer  in  the  limit 
0)  ->  0,  and  L/co  =  0(1)  is  shown  in  Fig.  3a  [2],  The  fuel-consumption  reaction  I  will 
occur  in  a  thin  layer  of  thickness  of  order  8co-2/3,  which  for  convenience  is  presumed  to  be 
located  at  £  =  x/8  =  0.  The  structure  of  this  layer  resembles  that  of  a  diffusion  flame  into 
which  fuel  diffuses  from  one  side  and  H-radicals  from  the  other.  Outside  this  diffusion 
layer  in  the  region  of  positive  x  there  exists  a  radical  nonequilibrium  layer  of  thickness  of 
order  Sw-1 ,  where  the  concentration  of  fuel  is  zero.  From  analyzing  the  structure  of  the 
diffusion-flame  layer  shown  in  Fig.  3a  it  can  be  shown  [2]  that  the  quantity  L  is  given  by 
the  expression 

L0  =  2co(l  +2k  +  40/3  +  8\j//15  +  2g),  (13) 

where 


K 


k5CM(k7  +  Yk8) 

kif(k6  +  k7+Yk8) 


o 


0  = 


YkibK  t  (Yk24  +  k2s)  Cm 
k2fXH2  klfX02 


YkibR 

f 

0 

,  o  = 

k2fk10bXH2  ^klbk18X02 
/ 

0 

,k2fXH2. 

n 

o 

X 

jT 

(m 

_ _ I 

, 

(14) 


A  schematic  illustration  of  the  presumed  structure  of  the  inner  layer  in  the  limit 
co  <*  with  L  =  0(1)  is  shown  in  Fig.  3b.  Here,  the  fuel  consumption  reaction  I  occurs 
in  a  relatively  broad  layer,  and  embedded  in  this  layer  is  a  thin  layer,  of  thickness  of  order 


5co  ^3,  where  reaction  IV  is  not  equilibrium.  Asymptotic  analysis  of  the  structure  of  this 
layer  shows  that  [2] 


-1  _  8 


L  1  = 


15 


Li* 


a 


2-  +  l|9  +  |v  +  |x  +  2a  +  UR 


3  -  4j32  +  p4 

p4 


ln(l  +  p) 


2p4-4p3-9p2  +  3P  +  6 


i-i 


2p3(l  +  P) 

The  various  parameters  appearing  in  Eq.  (15)  are  defined  as 

P  ■ 


-  n(ln|i  +  lnE^  +  C  -  l) 


Yk2bknR 

0 

1 1  — 

k2fkllbXH2 

0 

1  ,  kisXo2 

kifk^Xoj. 

,  - 

kifk13Xo2R_ 

1  T 

kllbXn2 

X  = 


kieRYpuW 

knWF 


=  (Ei  ib  +  nnbRTr  +  E2f+  n2fRTr-Eif)5(Tc-Tu)T0 

RT02 


(15) 


(16) 


where  En  and  nn  refer,  respectively,  to  the  activation  energy  and  the  temperature  exponent 
of  the  frequency  factor  of  the  elementary  chemical  reaction  n  shown  in  Table  1,  Tr  being  the 
reference  temperature,  set  equal  to  1600  K  [1,2],  and  C  =  0.5772  is  Euler's  constant. 

An  ad-hoc  approximation  to  determine  L  for  all  values  of  G)  has  been  proposed  and 
tested  previously  [2]  and  is  given  by  the  expression 


L  =  L„ 


1  + 


Lp 

0.18L.. 


(17) 


Equation  (17)  will  be  used  with  the  equations  derived  from  analyzing  the  structure  of  the 
H2-CO  oxidation  layer  to  calculate  the  burning  velocity  of  the  flame. 

6 .  Analysis  of  the  H2-CO  Oxidation  Layer 

The  burning  velocity  vu  can  be  calculated  from  Eqs.  (8)-(17)  if  T°,  X$2,  Xq2  and 
Xh2q  are  known.  To  determine  these  quantities  the  structure  of  the  H2-CO  oxidation  layer 
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downstream  from  the  inner  layer  must  be  analyzed.  In  this  layer  simplifications  to  the 
conservation  equations  arise  because  convection  can  be  neglected  in  the  first  approximation 
and  Xp  =  0;  hence  coi  =  0.  The  H  radicals  are  presumed  to  be  in  steady  state  in  this  layer, 
whence  from  Eq.  (5)  com  =  coiv-  It  can  then  be  shown  [2]  that  for  small  0  and  k,  Xh  is 
given  by  the  expression 

XH  =  R(l  -  0/2- k/2).  (18) 

The  water-gas  shift  reaction  II  is  presumed  to  be  in  partial  equilibrium  everywhere  in 
this  layer  except  in  a  thin  non-equilibrium  sublayer  adjacent  the  inner  layer.  Following  the 
development  in  Refs.  1  and  2,  the  structure  of  the  oxidation  layer  will  first  be  analyzed  by 
assuming  that  reaction  II  is  in  equilibrium  everywhere,  and  later  corrections  to  the  value  of 
X$2  and  Xq0  introduced  by  non-equilibrium  of  reaction  II  will  be  obtained.  Equilibrium  of 
reaction  II  implies  that  there  exists  an  algebraic  relation  between  the  concentrations  of  H2 
and  CO  that  can  be  written  as 

xco  =  ocxh2,  (19) 

where 

a  =  K3Xco2Lh2/(K9XH2oLco)  ,  xj  =  Xj/Li ,  (20) 


in  which  K9  =  0.0096  exp  (1 1623/T)  is  the  equilibrium  constant  of  elementary  reaction  9. 
The  partial-equilibrium  assumptions  make  it  necessary  to  eliminate  con  from  Eqs.  (5).  With 
com  =  0)iv  and  the  neglect  of  the  convective  terms  in  Eqs.  (5)  the  expressions 


[(1  +  a)xn2]  =  2o)ni  , 

dx2 


dL 

dx2 


[xh2  +  xhjq]  =  0  , 


77 [(1  +  «)xh2-  2xoJ  =  0  , 
dx2 


■fr  1(Qiii/2  +  a  Qrn/2  +  aQu)  xh2  +  't] 

dx2 


0, 


(21) 
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are  obtained,  where  Qni  =  Qra  +  Qiv-  Since  the  numerical  values  Qni/2  =  0.31 16  and 
Qnj/2  +  Qn  =  0.3479,  obtained  from  thermochemical  tables,  are  sufficiently  close,  they  can 
be  replaced  by  their  arithmetic  mean,  q  =  Qj^/2  »  Qjy/2  +  Qn  =  0.33.  The  expansions 

2qx  =  erj,  t  =  x^-ez,  q(l  +  a)xn2  =  e(b  +  z), 

2qxo2  =  e(a  +  z),  xj  =  xjb  -  ez;,  i  =  H2O,  CO2, 


are  then  introduced,  where  'tb  =  (Tb-Tu)/(Tc-Tu),  b  =  q(l  +  a)XH2b/(eLHz). 
a  s  2qXo2b/(eLo2)»  in  which  e  is  presumed  to  be  a  small  quantity,  and  T|,  z,  a,  b,  a,  z\t  Tb 
and  x  are  order  of  unity. 

Under  the  assumption  that  (yk24  +  k25)CM/(Yk'2b)  is  small,  an  expression  for  the  R  of 
Eq.  (11)  may  be  obtained  from  Eqs.  (3)  and  (22)  and  substituted  into  Eq.  (18)  to  show  that 
Xh  can  be  expressed  as 


XH 


e2K}/2Ky2K3lj^L1cg(z  +  b)3/2(z  +  a)1/2 
21/2q2XH2o(l  +  « f12 


9  _K 

2  2  J  ’ 


(23) 


where  Ki  =  12.7  exp  (-8108/T)  and  K2  =  2.25  exp  (-1040/T)  are  the  equilibrium 
constants  of  elementary  reactions  1  and  2,  respectively.  The  source  term  (Dm  in  Eq.  (21) 
can  be  written  from  Eqs.  (2)  and  (7)  in  terms  of  the  expansions  shown  in  Eqs.  (22)  as 

(Dm  =  2qe3Dm[Giii(z  +  a)3/2(z  +  b)3/2  +  GmS0(z  +  a)(z  +  b)2],  (24) 

where 


GO 
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S  = 


A(k5CM)Q(Kl/2Kl/2K,)0(LH,LoI)3/2(l  -  k/2  -  9/2) 
m  25/2q4XHjo(l  +  a°)3/2 

(2KiK2LH2)1/2(l  -  k/2  -  e/2)(k24  +  k25/Y) 

k5(l  +  a)1/2L^2 

G  (k5  Cm)  (k^  K  j 12  K2/2)  X&q  ( 1  +  a0)3/2(x/cp)  To 

111  (k5CM)0  (k3 K }/2K2/2)° XH2o ( 1  +  a)3/2(x./cp)°  T 

(KiK2)1/2(k24  +  k^/y)  kg(l  +  a°)1/2 
'(KiK2)1/2]°(k24  +  k25/Y)°k5  (1  +  a) 1/2  m 


Gm  = 


(25) 


In  the  inner  layer  Gffl  =  =  1,  and  following  previous  analysis  [1,2]  the  quantity  8  will 

be  presumed  to  be 

e  =  Dn1/4.  (26) 


Introduction  of  Eq.  (22)  into  Eq.  (21),  followed  by  use  of  Eqs.  (24)-(26)  with 
Gffl  =  GIU  =  1,  results  in  the  leading-order  problem 

=  (z  +  a)3/2(z  +  b)3/2  +  S°(z  +  a)(z  +  b)2,| 
dri2 

=  -1  at  T)  =  0, 


->  0  as  ri  ->  ■ 


which  implies  that 


)3/2(z  +  b)3/2  +  S°  (z  +  a)(z  +  b)2]dz  =  1 . 


(28) 


Here  a  =  a°  is  employed  in  the  definition  of  b  so  that,  like  a  and  S°,  this  parameter  may  be 
treated  as  constant  in  performing  the  integration.  It  can  be  easily  verified  that  if  a  =  b  =  0. 
then  this  problem  reduces  to  the  leading-order  problem  considered  earlier  [1,2]. 


ei 
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The  solution  to  Eq.  (28)  for  z°  as  a  function  of  a,  b  and  S°  must  be  obtained 
numerically.  For  a  given  value  of  <{>,  thermochemical  calculations  would  yield  the  values  of 
Tb,  Xo2b.  Xn2b,  Xcob>  Xco2b  and  Xn2ob-  The  quantities  a  and  b  depend  on  £,  which  may 
be  expressed  as  a  function  of  T°  according  to 


Tb-T° 
z°(Tc-Tu)  * 


(29) 


In  addition,  S°  depends  on  T°;  hence  all  parameters  in  Eq.  (28)  depend  on  results  from  the 
structure  of  the  inner  layer.  Before  calculating  the  burning  velocity,  corrections  to 
resulting  from  non-equilibrium  of  reaction  II  in  the  sublayer  adjacent  to  the  inner  layer  must 
be  obtained  [1,2].  In  this  layer  the  influence  of  reactions  III  and  IV  can  be  neglected. 
Introducing  the  expansions 

x  =  vC,  Xh2  =  Xh2  +  vLh2zh2,  Xco  =  Xco  +  vLcozco,  (30) 

where 

v  _  e^ksCM^Ka^ljf Ljf  (l  -  k/2  -  8/2)1/2 

2kffx§“  qwqg(zP  +  b)>'V  f  a)'«(l  +  a»)*'2  ’ 


it  can  be  shown  following  previous  [1,2]  methods  that  the  solution  in  the  non-equilibrium 
sublayer  can  be  written  as 


zh2  =  - 


1  -  a1 


(l  +  a°)3/2 


exp 


-{l  +a°)1/2^] — — 


1  +  a1 


(32) 


zco  = 


. . . exp  [  -  (l  +  a°)1/2C  ]  -  • 

(l  +  a°)3/2  1  +  a0 
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At  the  inner  layer  £  =  0;  hence 


e(z°  +  b) 
q(l  +  a0) 


vq(l  +  a°)  ~p  _  e(z°  +  b) 
e(z°  +  b)  H\  q(l+a°) 


1 


ci/2  (k5CMffl  -  a0)K30l\^2L^2(l  -  k/2  -  6/2)1/2 
Ikfxffil  +  a°}qi/2L^(zO  +  bP4(z  +  a)"4 


(33) 


which  can  be  used  in  the  expression 

L  ^l’fkiaLozll  +  a°)3/2(z0  +  a)5/2(x^q/E)3/2 
(k5  CM)° kj®  f  k't 3 LF ( 1  -  k/2  -  0/2) 


(34) 


obtained  with  the  aid  of  Eqs.  (8),  (10),  (11)  and  (23).  By  use  of  Eq.  (33),  Eq.  (34) 
expresses  L  in  terms  of  T°  and  z°,  since  e  is  related  to  T°  as  shown  in  Eq.  (29).  Since 
Eq.  (17)  provides  an  independent  expression  for  L  as  a  function  of  T°,  the  quantities  T° 
and  z°  can  both  be  calculated  numerically  when  Eq.  (28)  is  integrated.  In  view  of 
Eqs.  (26)  and  (29),  the  burning  velocity  can  then  be  calculated  by  rewriting  the  first  of 
Eqs.  (11)  and  (25)  as 


v.f  =  ■ 


Yfu 

A.]c 

>  fT 

Au 

WF 

LcpJ 

'j'0 

2  (k5  CM)°  !k1/2kV2K3)0(Lh,Lo?)3/2(i  -  K/2  -  9/2)  (Tb-  T°)4 
25/2q4Xj!jzo(l  +  a0)3/2(Tc-Tu)4z°4 


(35) 


7 .  Results  and  Discussions 

Equation  (35)  was  used  by  the  procedure  described  above  to  calculate  the  burning 
velocity  as  a  function  of  the  equivalence  ratio,  initial  temperature,  and  pressure.  In  these 
calculations  the  value  of  X./cp  appearing  in  Eq.  (35)  was  expressed  as  A,/cp  =  2.58  x 
10-4  (T/298)0-7  g/(cm  s).  The  Lewis  numbers  for  the  various  species  were  presumed  to 
be  constant,  with  Lp  =  0.97,  L02  =  1.1,  LhjO  =  0.83,  Loo*  =  1.39,  Lh2  =  0.3,  Lh  =  0.18 
and  Loo  =  1.11.  Calculations  were  performed  for  values  of  §  between  0.5  and  1.55,  for 
values  of  p  between  1  atm  to  40  atm,  and  for  Tu  =  300  K  and  Tu  =  600  K. 
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Figure  4  shows  results  of  numerical  solution  of  Eq.  (28)  for  the  value  of  z°  at  the 
inner  layer,  and  the  values  of  a  and  b.  The  rapid  decreases  of  a  on  the  rich  side  and  of  b  on 
the  lean  side  show  that  the  concentration  of  oxygen  in  the  post-flame  zone  is  small  for  rich 
flames,  and  the  concentration  of  hydrogen  in  the  post-flame  zone  is  small  for  lean  flames. 
Figure  5  shows,  the  values  of  L  and  co  as  functions  of  <)>,  demonstrating  that  co  decreases 
with  increasing  <)).  For  fuel-lean  flames,  where  the  value  of  to  is  large,  the  value  of  L  is 
given  approximately  by  Eq.  (15),  and  for  fuel-rich  flames,  where  the  value  of  co  is  small,  L 
is  given  approximately  by  Eq.  (13).  It  has  been  shown  previously  [2]  that  for 
stoichiometric  flames  the  value  of  L  decreases  and  the  value  of  co  increases  with  increasing 
pressure.  Hence,  the  expression  for  L  given  by  Eq.  (15)  becomes  more  accurate  for  fuel- 
lean  flames  at  high  pressure,  and  the  structure  of  such  flames  would  resemble  that  shown 
in  Fig.  3b.  Equation  (13)  would  become  more  accurate  for  fuel-rich  flames  at  low 
pressure,  and  the  structures  of  such  flames  resemble  that  shown  in  Fig.  3a.  However,  for 
fuel-rich  flames,  because  of  various  uncertainties  in  the  analysis,  such  as  questions 
concerning  our  neglect  of  C2  chemistry  in  the  chemical-kinetic  mechanism  shown  in 
Table  1,  the  extent  of  occurrence  of  the  diffusion-flame  layer  requires  further  study. 

In  Fig.  6,  the  variations  of  k,  0  and  a  with  <(>  are  shown,  and  in  Fig.  7  the  variations 
of  8,  v  and  e,  which  represent  the  thicknesses  of  the  various  layers  shown  in  Figs.  2,  3a 
and  3b,  are  plotted  as  functions  of  <j).  The  former  shows  the  current  expansions,  treating  k, 
0  and  a  as  small,  are  reasonably  accurate  except  for  rich  flames,  and  the  latter  shows  the 
relative  orderings  to  be  reasonable,  except  that  v  is  not  small  compared  with  £,  so  that 
water-gas  nonequilibrium  should  be  considered  throughout  the  £  layer  for  greater  accuracy. 
At  p  -  1  atm,  the  values  of  y  and  %  are  not  small  [2],  so  there  is  some  inaccuracy  in  the 
expansions  in  Eqs.  (13)  and  (15),  although  the  effects  tend  to  be  mitigated  by  other  terms; 
at  high  p  these  inaccuracies  disappear.  In  Fig.  8  the  temperature  in  the  inner  layer  T°  is 
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plotted  as  a  function  of  <J).  The  value  of  T°  increases  with  increasing  values  of  <{>  and 
approaches  Tb  for  very  fuel-rich  and  very  fuel-lean  flames.  Since  Eq.  (35)  shows  that  the 
burning  velocity  is  proportional  to  (Tb  -  T0)4,  the  burning  velocity  must  become  small  for 
highly  fuel-rich  or  fuel-lean  flames. 

In  Figs.  9  and  10  the  burning  velocity  vu  is  plotted  as  a  function  of  <{>  for  various 
values  of  p,  for  Tu  =  300  K  and  for  Tu  =  600  K,  respectively.  For  p  =  1  atm  and 
Tu  =  300K,  the  solid  curves  in  Fig.  9  indicate  that  the  burning  velocity  reaches  a 
maximum  value  of  approximately  36  cm/s  at  <j)  =  1.1.  Full  numerical  calculations  using 
somewhat  different  rate  data  [7]  gave  quite  similar  results,  except  that  the  maximum 
burning  velocity  was  much  closer  to  the  accepted  maximum  of  about  43  cm/s.  A  similar 
observation  applies  for  the  numerical  calculations  with  the  rate  data  of  Table  1 ,  as  shown  in 
Fig.  9.  In  addition,  it  is  seen  in  Fig.  9  that  for  p  =  1  atm,  <[>  =  1  and  Tu  =  300  K,  the  value 
vu  =  33  cm/s,  found  here  by  the  asymptotic  analysis,  lies  below  the  value  obtained  from 
detailed  numerical  calculations  employing  the  same  rate  parameters,  but  the  agreement 
improves  greatly  at  lean  conditions  and  at  higher  pressures.  As  the  pressure  increases,  the 
peak  value  of  vu  obtained  from  the  asymptotics  occurs  at  decreasing  values  of  $,  but 
always  greater  than  unity,  in  agreement  with  numerical  and  experimental  results. 

The  qualitative  features  of  the  behavior  of  vu  with  <j)  and  p  for  Tu  =  600  K  shown  in 
Fig.  10  are  similar  to  those  for  Tu  =  300  K  shown  in  Fig.  9,  but  the  values  of  vu  for  given 
values  of  <))  and  p  are  higher  for  Tu  =  600  as  is  expected.  It  is  interesting  to  note  that  for 
both  values  of  Tu  solutions  were  obtained  for  values  of  <[)  beyond  the  reported  [8,9]  rich 
and  lean  flammability  limits.  Therefore  the  results  of  the  analysis  indicate  that  there  does 
not  exist  a  chemical-kinetic  mechanism  for  flame  extinction  and  that  the  experimentally 
observed  flammability  limits  must  be  influenced  by  heat  losses. 

The  small  discontinuity  in  slope  at  <{>  =  1,  seen  in  Figs.  9  and  10  as  well  as  in  various 
other  curves,  is  a  consequence  of  the  expansion  about  Tc  in  the  analysis  of  the  layer  of  CO 
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and  H2  oxidation.  Figure  la  shows  the  discontinuity  in  the  slope  of  Tc  that  is  responsible 
for  this  behavior.  The  discontinuity  involves  terms  of  order  e  in  the  analysis  of  the 
oxidation  layer,  and  it  can  be  removed  by  alternative  treatments  of  this  layer.  For  example, 
expansions  can  be  pursued  about  Tb  instead  of  Tc,  with  xh2  (and  consequently  xco)  treated 
as  being  of  order  unity  in  Eq.  (22),  instead  of  being  of  order  e.  Because  of  the  magnitude 
of  e,  seen  in  Fig.  7,  these  different  approaches  give  results  differing  by  amounts  on  the 
order  of  20%,  as  has  been  verified  in  the  present  study  by  pursuit  of  alternative  analyses 
not  reported  here.  The  discontinuity  was  retained  in  the  present  paper  to  help  to  illustrate 
better  what  types  and  magnitudes  of  inaccuracies  may  be  anticipated  from  the  asymptotic 
treatment.  Future  research  could  be  directed  towards  seeking  optimum  approaches  to  the 
analysis  of  the  oxidation  layer,  with  higher-order  terms  in  e  considered  and  with  improved 
treatments  of  water-gas  nonequilibrium. 

It  has  been  shown  previously  [1,2]  that  an  effective  activation  energy  Eeff  may  be 


obtained  by  fitting  vu  to  an  Arrhenius  form  in  Tb-  Hence,  if  Eeff  =  2RTbd[ln(puvu)]/ dTb, 
then  from  differentiation  of  Eq.  (35)  with  respect  to  Tb  under  the  assumptions  that  T°  is 
independent  of  Tb  and  that  Tb  equals  Tc,  the  relationship 

T°-Tu 


Ecff  = 


4RTb 

Tb~Tu 


dTu 

|_Tb-T0  '  dTbJ 


(36) 


is  obtained.  If  the  quantity  dTu/dTb  is  neglected,  an  effective  Zel'dovich  number  [3]  can 
then  be  defined  as 

7  _  Eeff(Tb-Tu)2  _  4  ,07% 

RTj2(T°-Tu)  ezo  ’ 


where  use  was  made  of  Eq.  (29).  In  Fig.  11,  this  Ze  is  plotted  as  a  function  of  <j)  for 
various  values  of  p  at  Tu  =  300  K.  Interestingly,  the  value  of  Ze  is  a  minimum  at  4>  =  1.0, 
and  increases  rapidly  with  increasing  or  decreasing  values  of  <)).  In  addition,  the  minimum 
value  of  the  Ze  increases  with  increasing  values  of  p.  The  large  value  of  Ze  near  the 
experimentally  observed  flammability  limits  implies  that  the  flame  is  extremely  sensitive  to 
heat  losses,  and  hence  it  would  be  increasingly  difficult  to  obtain  a  steady  flame 
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propagation  in  increasingly  fuel-lean  or  fuel-rich  mixtures.  These  observations  could  bear 
on  why  flammability  limits  are  observed  in  experiments.  It  must  be  emphasized,  as 
discussed  in  Ref.  1,  that  the  large  value  of  the  Zel'dovich  number  obtained  here  is  related  to 
T°,  which  is  determined  by  the  relative  rates  of  important  elementary  reactions,  and  is 
unrelated  to  one-step  activation-energy  asymptotics. 

8.  Conclusions 

This  paper  extends  the  asymptotic  analysis  of  methane  flames  initiated  in  [1]  and  [2] 
by  considering  off-stoichiometric  conditions.  It  shows  that  such  an  analysis  in  principle 
can  be  done  but  that  many  parameters  enter  into  the  formulation.  Nevertheless,  the  essence 
of  the  structure  originally  proposed  in  [1]  remains  valid  in  the  entire  range  of  equivalence 
ratios  and  pressures  considered  here.  A  particularly  useful  result  is  the  possibility  of 
defining  and  calculating  an  effective  Zel'dovich  number  and  thereby  establishing  a  link  to 
previous  large-activation-energy  analyses. 
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Rate  constants  are  kn  =  AnTn~exp(—HJn/\RT))]  units  are  moies,  cubic  centimeters,  sec¬ 
onds,  degrees  Kelvin,  and  calories/rnole. 

a)  Third  body  efficiencies  are:  CII^  :  C.5,  II20  :  6.5,  C02  :  1.5,  II2  :  1.0,  CO  :  0.75,  02  : 
0.4,  N2  :  0.4,  all  other  species:  1.0. 

b)  High  pressure  vaiue  Th.e  pressure  dependence  is  given  by  the  Lindemann  form 
k  =  A.”oo/(l  4-  (oilT/;>)),  where  aR  =  0.517exp(— 9000/T).  with  p  in  atm  and  T  In 
degrees  Kelvin. 
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CHAPTER  VI. 


THE  STRUCTURE  AND  EXTINCTION  OF  PARTIALLY  PREMIXED  FLAMES 

BURNING  METHANE  IN  AIR 

It  was  mentioned  in  chapter  II  that  turbulent  reacting  flows  can  be 
modelled  as  a  statistical  ensemble  of  laminar  diffusion  flamelcts,  laminar 

premixed  flamelets,  and  laminar  partially  premixed  flamelets.  These  flamelets 
are  subjected  to  strain.  In  chapters  III  and  IV  the  structure  of  diffusion 

flamelets  were  considered,  while  in  chapter  V  the  structure  of  premixed 
flamelets  were  considered.  To  complete  the  description  of  laminar  flamelets, 
in  this  chapter  laminar  partially  premixed  flamelets  are  considered.  As  in 
chapters  III,  IV  and  V,  methane  is  used  as  the  model  hydrocarbon  fuel.  A 
experimental  and  numerical  study  of  the  structure  and  mechanisms  of 

extinction  of  counterflow,  partially  premixed,  flames  burning  methane  was 

performed.  Experimental  data  was  obtained  for  the  temperature  profile,  the 
axial  velocity  profile,  the  concentration  profiles  for  various  stable  species  and 

the  critical  conditions  of  the  flame  at  extinction.  The  results  WErc  compared 
with  numerical  calculations  performed  at  conditions  identical  to  those  used  in 
the  experiment.  Experimental  measurements  and  numerical  calculations  show 
that  the  partial  prcmixi'ng  of  the  reactant  streams  of  a  counterflow  diffusion 

flame  makes  the  flame  less  resistant  to  stretch.  In  addition,  numerical 
calculations  show  that  the  reaction  zone  of  a  partially  premixed  flame  exhibits 
a  diffusion  flame  and  a  premixed  flame-like  structure.  Both  results  are 
consistent  with  previous  asymptotic  analyses  of  the  structure  of  partially 

premixed  flames. 

The  research  summarized  above  has  been  published  in  the  Proceedings  of 

the  Twenty-Second  Symposium  (International)  on  Combustion,  The  Combustion 
Institute,  pp  1555-1563,  1988.  The  coauthors  of  this  manuscript  arc  Professor 
M.  D.  Smooke,  and  Dr.  I.  K.  Puri.  The  research  performed  by  Professor  M.  D. 
Smookc  at  Yale  University  was  supported  by  the  Office  of  Naval  Research. 


CHAPTER  VII.  HYDROGEN-AIR  DIFFUSION  FLAMES 


VII.l  ANALYSIS  OF  THE  STRUCTURE  OF  COUNTERFLOW  HYDROGEN- 

AIR  DIFFUSION  FLAMES. 

Numerical  calculations  and  asymptotic  analysis  of  the  structure  of 
hydrocarbon  flames  described  in  chapters  II  -  VI  show  that  the  structure  of 
the  oxidation  layer  in  these  flames  are  primarily  determined  by  H2-O2 
reactions.  Hence,  to  obtain  an  improved  understanding  of  the  structure  of  the 

oxidation  layer  in  hydrocarbon  flames,  numerical  calculations  were 

performed  to  determine  the  structure  of  diffusion  flames  stabilized  between 
counterflowing,  laminar  streams  of  hydrogen  and  air.  The  calculations  were 
performed  at  different  values  of  the  rate  of  strain  and  at  a  value  of  the 
absolute  pressure  equal  to  one  atmosphere.  The  detailed  kinetics  of  oxidation 
of  hydrogen  was  represented  by  eleven  elementary  reactions  involving  seven 

species.  Temperature  profiles,  concentration  profiles  of  various  stable  species 
and  radicals,  rates  of  various  elementary  chemical  reactions,  and  rates  of 
formation  of  various  species  and  radicals  were  calculated  for  a  number  of 
values  of  the  rate  of  strain.  Results  are  plotted  at  a  low  value  of  the  rate  of 

strain  (=  30  s'1)  and  at  a  high  value  of  the  rate  of  strain  (=  1440  s'1).  At  the 

low  value  of  the  rate  of  strain,  a  number  of  reactions  were  found  to  be  in 

partial  equilibrium  over  a  major  part  of  the  flow  field.  Using  approximations 
suggested  by  the  results  of  the  numerical  calculations,  the  equilibrium 

structure  of  the  flame  was  determined  and  the  results  were  compared  with  the 
results  of  numerical  calculations  performed  at  a  value  of  the  rate  of  strain 
equal  to  =  30  s'1  and  good  agreement  was  obtained  for  the  concentration 
profiles  of  H2,  O2,  H2O,  and  H.  However  the  shapes  of  the  temperature  profiles 
were  different  and  they  are  attributed  to  nonequilibrium  of  the  major  chain 
branching  reaction  H  +  O2  =  OH  +  0.  This  study  shows  that  a  number  of 
approximations  employed  in  chapters  I  -  VI  in  describing  the  structure  of  the 

oxidation  layer  of  hydrocarbon  air  flames  arc  reasonably  accurate. 

The  research  summarized  above  has  been  accepted  for  publication  in 
Progress  in  Astronautics  and  Aeronautics,  1990.  The  coauthors  of  t  he 
publication  arc  Dr.  Vcnkat  Tangirala,  Professor  C.  Trevino,  and  Professor  M.  D. 
Smooke. 


VII.2  THE  INFLUENCE  OF  THE  LEWIS  NUMBER  OF  THE 
REACTANTS  ON  THE  ASYMPTOTIC  STRUCTURE  OF  COUNTERFLOW 
AND  STAGNANT  DIFFUSION  FLAMES 


In  the  asymptotic  analysis  described  in  chapter  VII.  1,  the  Lewis  number 
for  all  species  were  presumed  to  be  unity.  However,  for  hydrogen  the 
characteristic  Lewis  number  is  small.  Hence,  in  this  section  the  asymptotic 
structure  of  counterflow  and  stagnant  diffusion  flames  are  analyzed  in  the 
limit  for  large  values  of  the  overall,  nondimensional  activation  energy,  Ta, 
characterizing  the  rate  of  the  reaction,  and  results  are  given  for  small  values 
of  the  stoichiometric  fuel  to  oxygen  mass  ratio.  The  chemical  reaction 
between  the  fuel  and  the  oxidizer  is  represented  by  a  one-step,  irreversible 
process.  A  new  approach  is  developed  to  characterize  the  influence  of  the 
Lewis  number  of  the  fuel,  Lf,  and  the  Lewis  number  of  the  oxidizer,  Lo,  on  the 
outer  and  the  inner  structure  of  near  equilibrium  diffusion  flames.  Explicit 
algebraic  formulas  to  predict  the  critical  conditions  of  flame  extinction  are 
also  given. 

For  counterflow  diffusion  flames  at  fixed  values  of  Lo,  the  flame  moves 
significantly  toward  the  oxidizer  stream,  and  the  heat  losses  toward  the 
oxidizer  region  of  the  flame  increases  significantly  with  decreasing  values  of 
Lf-  The  value  of  the  maximum  flame  temperature  is  relatively  insensitive  to 
the  variations  in  Lf  although  the  value  of  the  rate  of  strain  at  extinction,  A, 
increases  significantly  with  decreasing  values  Lf  and  increasing  values  of  Ta. 
At  fixed  values  of  Lf  and  decreasing  values  of  Lo ,  the  flame  moves  slightly 
toward  the  fuel  stream;  the  heat  losses  toward  the  fuel  stream  increases 
slightly,  and  there  is  moderate  increase  in  the  value  of  the  maximum  flame 
temperature.  The  value  of  A  increases  with  decreasing  values  of  Lo  for  large 
values  of  Ta  and  is  relatively  insensitive  to  variations  in  Lo  for  moderate 
values  of  Ta. 

The  inner  and  outer  structure  for  stagnant  diffusion  flames  where 
convection  is  absent  arc  qualitatively  similar  to  those  for  counterflow 
diffusion  flames.  However,  the  value  of  the  maximum  flame  temperature 
increases  significantly  with  decreasing  values  of  Lq  and  fixed  values  of  Lf. 


The  results  developed  here  are  used  to  obtain  overall  chemical  kinetic 
rate  parameters  characterizing  the  gas  phase  oxidation  of  methane  using 
previously  measured  values  of  the  critical  conditions  of  flame  extinction. 

The  research  summarized  above  has  been  published  in  Combustion 
Science  and  Technology,  Volume  64,  pp  243-261,  1989  The  coauthor  of  this 
publication  is  Professor  C.  Trevino. 


